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SUMMARY 


The  evaluation  of  any  propeller -driven  V/STOL  configuration  requires  analytical 
methods  that  adequately  assess  the  aerodynamic  characteristics  and  the  forces 
arising  from  the  propeller  and  the  slipstream. 

This  study  was  to  develop  a  unified  analytical  procedure  to  evaluate  the  effects  of 
passive  high-lift  devices  on  deflected-slipstream  or  tilt -wing  V/STOL  configura¬ 
tions.  Methods  were  developed  to  predict  the  two-dimensional  flapped  airfoil  char¬ 
acteristics  to  be  used  in  a  span  load  program.  The  span  load  results  are  used  in 
procedures  for  estimating  the  coefficients  of  lift,  longitudinal  force,  and  moment 
for  a  wing  partially  immersed  in  a  propeller  slipstream.  These  characteristics 
can  then  be  used  in  a  performance  program  developed  to  calculate  the  takeoff,  land¬ 
ing  and  transition  maneuvers.  In  addition  to  these  tasks,  investigations  were  made 
into  downwash  characteristics,  wind  tunnel  wall  corrections,  and  correlations  of 
flight  test  data  with  theory.  An  analysis  of  the  effects  of  high-lift  devices  on  the 
.  performance  of  a  tilt -wing  V/STOL  configuration  is  included  in  the  appendix. 

The  procedures  to  predict  the  two-dimensional  aerodynamic  characteristics  of  high- 
lift  devices  furnish  reasonable  estimates,  generally  within  10  percent  of  experi¬ 
mental  results,  that  can  be  used  in  the  span  load  program.  The  trend  from  General 
Dynamics  experimental  data  indicates  that  the  maximum  lift  levels  from  NACA  data 
are  approximately  10  to  20  percent  lower.  This  may  be  attributed  to  the  General 
Dynamics  testing  technique  which  utilizes  side-wall  blowing  to  eliminate  adverse 
wall  boundary  layer  interference  effects. 

The  program  initiated  at  Convair  to  develop  a  more  basic  approach  to  the  problem 
of  estimating  the  lift  and  longitudinal  force  coefficients  was  well  justified.  It  is 
considerably  less  limited  than  existing  empirical  methods  in  that  it  requires  only 
two-dimensional  data  for  application  to  general  configurations.  The  procedures  for 
predicting  lift  and  longitudinal  force  coefficients  of  a  wing-flap  combination  give 
satisfactory  results  at  all  thrust  coefficients.  However,  the  procedures  for  pre¬ 
dicting  pitching  moment  coefficients  result  in  erratic  correlations  for  cases  with 
flaps  deflected. 

The  entire  task  of  correlating  wind  tunnel  data  that  has  been  corrected  for  wall 
effects  with  actual  flight  test  conditions  could  not  be  accomplished  due  to  the  lack 
of  a  good  set  of  flight  test  data  for  comparison.  The  evaluation  of  tunnel  wall 
effects  indicated  that  currently  available  correction  procedures  for  lift  and  drag 
yield  erroneous  results  and  should  not  be  applied. 


ili 


The  complexity  of  the  performance  program  made  numerous  assumptions  and  capa¬ 
bility  restrictions  desirable.  For  example,  the  equations  of  motion  were  reduced 
to  two  dimensions  (i.e. ,  no  pitch  dynamics),  and  several  terms,  such  as  pitch 
rates,  wing  tilt  angular  rates  and  pilot  braking  response  time,  were  used  as  pre¬ 
determined  inputs.  These  assumptions  were  made  to  expedite  the  development  of 
the  system;  however,  the  accuracy  of  the  methods  employed  is  sufficient  to  reflect 
the  impact  of  high-lift  devices  on  V/STOL  aircraft.  The  program  is  a  viable  tool 
for  performance  estimation. 

The  problems  associated  with  predicting  the  aerodynamic  forces  acting  on  a  wing 
immersed  in  a  slipstream  should  be  investigated  further.  To  improve  the  basic 
procedures,  it  is  necessary  that  (1)  additional  correlations  of  pitching  moments  be 
made  with  experimental  data  to  develop  improved  empirical  factors  for  large  chord 
ratio  flaps,  (2)  methods  be  developed  to  incorporate  the  pronounced  nonlinear  char¬ 
acteristics  of  leading  edge  devices  at  low  angles  of  attack,  (3)  the  downwash  esti¬ 
mating  procedure  be  programmed  for  digital  computer  application  and  correlations 
be  made  with  experimental  data,  and  (4)  tilt-prop  configuration  methodology  be 
developed  and  incorporated  into  the  aerodynamic  and  performance  programs. 

Jet  flap  and  boundary  layer  control  devices  should  be  investigated  so  that  methods 
could  be  developed  to  describe  the  effects  of  these  devices  on  a  wing  immersed  in 
a  slipstream. 

A  study  should  be  performed  to  define  a  wind  tunnel/flight  test  program  to  corre¬ 
late  V/STOL  vehicle  aerodynamic  characteristics. 
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FOREWORD 


This  report  presents  the  results  of  an  investigation  of  the  effects  of  high-lift  de¬ 
vices  on  V/STOL  aircraft  performance  and  stability  and  control. 
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1.0  INTRODUCTION 


Interest  in  V/STOL  aircraft  has  emphasized  the  need  to  improve  the  methods  of 
predicting  the  aerodynamic  characteristics  of  these  vehicles  at  very  low  forward 
flight  velocities.  A  general  understanding  of  the  high-lift  characteristics  of  cer¬ 
tain  configurations  exists  because  of  the  amount  of  available  experimental  data. 
However,  reliable  analytical  methods  of  predicting  the  high -lift  characteristics 
over  a  full  range  of  variables,  i.e. ,  angle  of  attack,  flap  angle,  thrust  coefficient, 
etc. ,  are  not  available. 

A  comprehensive  literature  search  on  V/STOL  high-lift  devices  was  conducted 
and  is  documented  in  Volume  II  of  this  report.  The  avaiUble  two-dimensional 
experimental  data  were  reviewed,  and  selected  NACA  test  data  and  data  developed 
by  General  Dynamics  were  used  in  developing  a  procedure  to  predict  the  two- 
dimensional  characteristics  of  the  following  passive  high -lift  devices: 

Leading  Edge  Devices 

Nose  Flaps 
Slats 
Kruegers 

These  characteristics,  along  with  experimental  clean  airfoil  data  from  Reference 
1,  are  used  in  a  span  load  program  to  predict  the  three-dimensional  characteris¬ 
tics  of  an  arbitrary  configuration. 

The  limitations  imposed  by  the  available  semiempirical  methods  for  estimating 
the  forces  on  a  wing  immersed  in  a  propeller  slipstream  indicated  the  require¬ 
ment  for  a  more  basic  approach  to  the  problem.  The  approach  developed  at 
Convair  depends  entirely  on  the  availability  of  two-dimensional  data  and  was  first 
conceived  by  Canad^r  Limited  in  the  early  1960's.  Thrust,  free-stream,  and  slip¬ 
stream  effects  on  the  unpowered  lift,  drag,  and  pitching  moment  are  predicted  by 
the  theory  discussed  in  this  report. 

These  data  are  then  trimmed  and  used  in  the  performance  program  to  integrate 
two-degree-of-freedom  equations  of  motion  to  predict  vertical  or  rolling  takeoff, 
landing,  and  transition  maneuvers. 

This  report  describes  the  s>tudy  effort  and  the  development  of  a  methodology  to 
evaluate  the  effects  of  high-lift  devices  on  V/STOL  aircraft  performance. 


Trailing  Edge  Flaps 


Plain  Flaps 
Single -Slotted  Flaps 
Double -Slotted  Flaps 
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TWO-DIMENSIONAL  AERODYNAMIC  CHARACTERISTICS 
OF  HIGH-LIFT  DEVICES 


The  design  of  V/STOL  aircraft  requires  procedures  that  adequately  predict  the 
aerodynamic  characteristics  of  high-lift  devices.  These  devices  are  usually 
movable  portions  of  the  leading  or  trailing  edge  of  the  wing  that  are  deflected  to 
increase  the  maximum  lift  coefficient. 

Considerable  data  are  available  from  NACA  sources  on  two-dimensional  tests 
that  were  conducted  during  the  1930’s.  Unfortunately,  these  tests  were  not 
planned  to  investigate  systematic  parameters  on  high-lift  airfoils.  Several 
attempts  have  been  made,  such  as  References  2,  3,  and  4,  to  organize  the  avail¬ 
able  data  in  a  useable  form.  The  above  references  either  treated  the  increment 
in  lift  below  stall  or  attempted  to  handle  the  complete  lift  curve. 

This  section  summarizes  a  procedure  (based  on  thin  airfoil  theory  where  possible) 
from  Reference  5  that  predicts  the  two-dimensional  aerodynamic  characteristics 
of  the  following  high -lift  devices: 

Trailing  Edge  Flaps  Leading  Edge  Devices 

Plain  Flaps  Nose  Flaps 

Single -Slotted  Flaps  Kruegers 

Double -Slotted  Flaps  Slats 

Selected  NACA  test  data,  mainlv  the  23  series  flapped  airfoils,  and  data  devel¬ 
oped  by  General  Dynamics  were  utilized  in  developing  the  procedure  to  predict 
incremental  effects  of  high -lift  devices.  The  predicted  increments  are  added  to 
the  clean  airfoil  characteristics  determined  either  from  the  experimental  data  in 
Reference  1  or  predicted  with  the  method  of  Reference  2. 

2.1  TRAILING  EDGE  FLAPS 

2.1.1  Theory 

Descriptions  of  the  basic  theoretical  treatment  of  the  effects  of  flaps  on  the  char¬ 
acteristics  of  airfoils,  using  an  extension  to  thin  airfoil  theory,  are  given  in 
References  6,  7,  and  8.  These  analyses  lead  to  expressions  by  which  lift  and 
pitching  moment  could  be  calculated.  A  further  refinement  of  the  theory,  in 
References  9  and  10,  shows  equations  relating  maximum  lift  increment  to  lift 
increment  at  a  =0.  These  equations  are  based  on  the  assumption  that  leading 
edge  separation  is  dependent  only  on  the  additional  "angle  of  attack"  loading. 
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Only  a  cursory  description  of  thin  airfoil  theory  is  included  since  it  is  adequately 
covered  in  a  number  of  references.  Consider  the  flapped  airfoil  shown  below: 


Flapped  Airfoil  Geometry 

The  angle  of  attack  of  the  chord  line  with  the  flap  deflected  is 


“eff 


-  -.,0 


From  thin  airfoil  theory,  the  lift  coefficient  is 

c .  =  2tr  A  +  ttA 
l  o  1 

where 


and 


ir-6f 

A  =  a  + - 5 

o  ff  f 


26f 

A  =  - sin  0 

1  IT  f 


The  equation  for  the  lift  coefficient  of  a  flapped  airfoil  may  be  written  as 

/  *"®f  \ 

CA  =  2ir  (a  +  — ~  6{)  +  26{  sinQf 


(1) 


(2) 


(3) 


(4) 
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cA  =  2ir a  +  2  (it -6{  +  sin0f)  6f 


The  two-dimensional  lift  coefficient  can  then  be  written  as 


=  2ira  +c^  6f  =  2rt(a-a^6i) 
6 


where 


cl  =  2  (tt-  0f  +  sin  0^) 


a6  ~  217 


The  theoretical  lift  effectiveness,  c^  ,  is  given  in  Figure  1.  The  rate  of  change 

of  zero-lift  angle  with  flap  deflection,  a *,  may  readily  be  obtained  from  c «  . 

0  6 

In  thin  airfoil  theory  the  ideal  angle  of  attack,  o^,  can  be  defined  as  the  angle  of 
attack  where  lift  is  obtained  from  camber  alone  with  no  suction  peak  at  the  nose . 
This  definition  results  by  considering  the  aerodynamic  loading  over  a  thin  cam¬ 
bered  airfoil  in  two  parts.  One  part,  the  basic  load  distribution,  is  characteristic 
of  the  camber -line  shape.  The  other  part,  the  additional  load  distribution,  is  due 
only  to  angle  of  attack. 


The  suction  peak  then  depends  on  angle  of  attack,  and  for  the  symmetrical  air¬ 
foil,  oqc  =  0.  The  symmetrical  airfoil  stalls  at  ag  and  the  cambered  airfoil 


would  stall  at  a  ,  where 
sc 


as  =  “sc“«ic 


It  follows  that  for  the  flapped  airfoil,  an  ideal  angle  of  attack  may  be  defined  such 
that 


a  13  a  -a., 
s  sf  if 


The  ideal  angle  of  attack  for  the  flapped  airfoil  must  be  such  that  the  AQ  in  Equa¬ 
tion  (3)  goes  to  zero.  It  follows  that  for  the  flapped  airfoil, 
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Figure 


2.  Theoretical  Values  of  the  Ratio  of  the  Increment  of  Maximun 
to  the  Increment  of  Lift  at  a=0  Versus  Flap 
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Substituting  in  Equation  (13),  we  may  write 


Ac  =  -  1/2  (sin  0  -  1/2  sin  20  ) 
m  f  f  I 


=  - 1/2  (sin  0f  -  sin  0f  cos  0f)  6{ 


where 


C°s0f  =  -1-2—) 


The  increment  of  moment  due  to  flap  deflection  can  now  be  expressed  as 


115) 


Ac  =  -  (l  -  sin  0  6, 
m  \  c  /  f  f 


f 


sin  0  6  =  c  6, 
c  f  f  m  f 
6 


(16) 


where  c  =  the  moment  effectiveness  parameter, 
m6 

For  convenience  in  analysis,  the  pitching  moment  increment  caused  by  flap  deflec¬ 
tion  can  be  expressed  as  the  following  ratio: 


Ac 

m 


Ac 


(17) 


The  theoretical  values  for  cm^/c^ 


are  given  in  Figure  3. 


2.1.2  Lift  Increment  at  a  =  0  Degrees 


The  lift  effectiveness  of  simple  trailing  edge  flaps  can  be  defined  from  thin  airfoil 
theory.  The  rate  of  change  of  lift  with  flap  deflection  at  constant  angle  of  attack  as 
given  by  Equation  (7)  is 
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The  lift  of  the  flapped  airfoil  is  related  directly  to  the  effective  airfoil  angle  of  attack 
for  the  linear  lift  range. 

The  theoretical  lift  effectiveness  (c^  )  is  primarily  a  function  of  flap  chord  and  air¬ 
foil  thickness  ratio.  Any  increase  in  airfoil  thickness  should  increase  the  theoreti¬ 
cal  lift  effectiveness  term.  However,  it  has  been  shown  experimentally  that  the 
boundary  layers  on  thick  airfoils  have  a  larger  effect  than  on  thin  airfoils  and  cause 
a  reduction  in  lift  effectiveness.  This  reduction  in  lift  effectiveness  due  to  viscous 
effects  is  greater  than  the  theoretical  increases  due  to  thickness. 

Since  the  major  effects  are  from  the  boundary  layer,  an  empirical  correlation 
parameter  has  been  developed  which  relates  the  lifting  characteristics  of  various 
trailing  edge  flaps  to  the  theoretical  lifting  values  for  a  flat  plate. 

This  parameter,  the  effective  turning  angle,  is  defined  by  the  following  equation: 


where 


(18) 


Airfoil  Trailing  Edge 

The  effective  turning  angle  is  used  to  relate  the  lifting  efficiency  of  the  flap  to  the 
theoretical  thin  airfoil  lift  values.  The  lifting  efficiency  is  determined  for  plain, 
single-  and  double -slotted  flaps  shown  in  Figure  4. 

The  procedure  assumes  no  compressibility  or  Reynolds  number  effects.  The  Rey¬ 
nolds  number  effect  is  accounted  for  in  predicting  the  characteristics  of  the  clean 
airfoil.  Reference  11  indicates  that  the  lift  increment  due  to  flaps  is  essentially 
constant  with  Reynolds  number.  The  data  used  in  developing  the  procedure  was 
generally  taken  from  two-dimensional  tests  conducted  at  Reynolds  numbers  of  the 
order  of  3.0  x  10 6. 
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The  lift  increment  at  a  =  0  can  then  be  calculated  for  plain  flaps  by  using  the  fol¬ 
lowing  expression: 

Ac*  =  np  ol  6f  (19) 

a  =  0  6 

where  r)p  is  the  lifting  efficiency  of  the  plain  flap  developed  from  data  in  Refer¬ 
ences  12  through  17  (from  Figure  5). 

Cjj  is  the  theoretical  lifting  effectiveness  (from  Figure  1). 

6 

is  the  flap  deflection  angle . 

Comparisons  of  low-speed  test  data  with  predicted  Ac^  show  that  the  corre¬ 
lations  generally  were  within  plus  or  minus  10  percent  of  the  experiment  values.  * 

The  introduction  of  a  slot  to  the  plain  flap  creates  slot  flow  that  energizes  the  bound¬ 
ary  layer  over  the  flap.  This  flow  at  the  slot  lip  decreases  the  viscous  effects 
present  on  the  plain  flap  and  therefore  increases  the  lifting  efficiency  of  the  single- 
slotted  flap  for  deflections  greater  than  15  degrees.  The  lift  increment  at  a  =  0  due 
to  flap  deflection  for  the  single-slotted  flap  is  given  by  the  following  expression: 

Act  =  \  °t  6f  (20> 

a=0  6 

Since  most  slotted  flaps  extend  the  airfoil  chord,  the  equation  is  further  modified 
as  follows: 


Ac 


j l 

a~0 


(21) 


where  r)^  is  the  lifting  efficiency  of  the  single -slotted  flap  developed  from  data  in 
References  18  through  27  (from  Figure  5). 

c  i  is  the  theoretical  lifting  effectiveness  based  on  flap  chord  (from 
61 

Figure  1). 

6f  is  th  :  single -slotted  flap  deflection. 

is  the  ratio  of  extended  chord  to  airfoil  chord  (from  geometry  in 
Figure  4). 

♦Note:  All  data  correlations  are  shown  in  Reference  5. 
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All  of  the  flaps  analyzed  incorporated  chord  extension  which  varied  from  0  to  ap¬ 
proximately  10  percent.  The  procedure  described  above  showed  very  good  results 
for  the  single -slotted  flaps  at  any  arbitary  deflection  and  chord  extension.  In  addi¬ 
tion,  single -slotted  flaps  with  100  percent  flap  chord  extension  (chord  extension 
equal  to  flap  chord)  were  analyzed. 

Because  of  the  lack  of  systematic  NACA  experimental  data  for  the  various  geometric 
and  aerodynamic  variables  involved,  it  is  very  difficult  to  arrive  at  a  more  accurate 
data  correlation  basis.  The  data  correlation  for  single -slotted  flaps  relied  heavily 
on  systematic  comparisons  of  data  from  References  18  through  27.  The  correla¬ 
tions  were  generally  within  plus  or  minus  10  percent  of  the  experimental  values. 


The  introduction  of  a  secondary  slot  to  the  flaps  is  handled  by  adding  a  term  to 
Equation  (21).  The  lift  increment  at  a  =0  due  to  flap  deflection  for  double-slotted 
flaps  is  then  given  by  the  following  expression: 


V°i 


Ac*  „  =  Vi.  6f,  \—  +  Vi  \  l1  +  —r-  <22 

a=0  6  1  \  /  60  2  \  / 

J-  u 

where  subscripts  1  and  2  refer  to  the  primary  and  secondary  flap,  respectively. 


\  c  /  is  the  ratio  of  extended  chord  to  airfoil  chord  resulting  from  deflec¬ 
tion  of  the  primary  flap. 


V°iV 


is  the  ratio  of  extended  chord  to  airfoil  chord  resulting  from 


deflection  of  the  secondary  flap. 


The  flap  geometry  is  from  Figure  4. 

The  secondary  slot  caused  by  deflecting  the  aft  flap  is  not  as  effective  in  creating 
lift  as  the  primary  slot.  As  a  consequence,  an  additional  empirical  correlation 
parameter  rj2 g  was  developed.  Values  of  versus  aft  flap  deflection  are  shown 
in  Figure  6.  As  the  forward  flap  deflection  is  increased  to  approximately  20  de¬ 
grees,  the  efficiency  of  the  aft  flap  decreases  for  aft  flap  deflection  greater  than 
15  degrees.  The  lifting  effectiveness  for  the  aft  flap,  rj^,  is  defined  as 


\  \  TJ2S 


(23) 


where 


Tj  is  determined  from  the  aft  flap  deflection  and  chord  (from  Figure  5). 
17  accounts  lor  reduced  effectiveness  of  the  aft  flap  (from  Figure  6). 
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A  similar  procedure  was  also  developed  to  analyze  the  double -slotted  flap  that  has 
a  vane  and  aft  flap.  The  vane  is  contoured  to  produce  slot  flow  at  the  forward  slot 
lip  and  at  a  forward  secondary  slot.  This  type  of  double -slotted  flap  is  shown  in 
Figure  4.  The  available  data  indicated  that  the  aft  flap  for  this  configuration  does 
operate  at  full  efficiency.  The  lift  increment  at  a  =  0  due  to  flap  deflection  for  a 
vane  and  aft  flap  combination  is  given  below: 


+  T)  c 


1  i  uf 


(24) 


where  subscripts  v  and  2  refer  to  the  vane  and  secondary  flap,  respectively, 
is  the  turning  efficiency  presented  in  Figure  5. 


<c  +c 


is  the  ratio  of  extended  chord  to  airfoil  chord  resulting  from 


deflection  of  the  vane. 


f  c  -c  \ 

2  v| 

1  + - 1  is  the  ratio  of  extended  chord  to  airfoil  chord  resulting  from 


deflection  of  the  secondary  flap. 


The  comparisons  of  experimental  data  from  References  18,  24,  and  27  through  33 
with  predicted  values  for  the  double-slotted  flaps  indicated  correlations  within  plus 
or  minus  10  percent  of  the  experimental  values. 


2.1,3  Lift  Curve  Slope  Increment 

The  theories  for  thin  airfoils  indicate  that  the  lift  curve  slopes  of  a  flat  or  cambered 
airfoil  are  the  same,  c^  =217.  The  camber  term  simply  shifts  the  lift  curve  upward 

ry  "*  “ 

by  a  constant  increment.  The  experimental  lift  curve  slopes  for  the  unflapped  and 
flapped  airfoils  without  chord  extension  generally  verify  the  theory.  This  assump¬ 
tion  is  valid  up  to  the  point  where  the  flow  separates  from  the  airfoil  or  flap.  After 
the  flow  separates,  the  lift  curve  slope  of  the  flapped  airfoil  decreases  to  a  value 
below  that  of  the  basic  airfoil. 


The  effect  of  flap  extension  is  determined  by  simply  considering  the  change  in  effec¬ 
tive  chord.  For  a  flapped  airfoil  with  chord  extension,  the  lift  curve  slope  is  then 
directly  related  to  the  airfoil  chord  increase. 
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c, 


(25) 


“EXT 


a 


AIRFOIL 


This  procedure  assumes  that  (1)  the  lift  curve  slope  for  plain  flaps  is  that  of  the 
unflapped  airfoil,  and  (2)  the  change  in  lift  curve  slope  for  flaps  that  extend  the  air¬ 
foil  chord  is  given  by 

=  c, 
l  l 

“ext  “airfoil 

where  c«  is  the  lift  curve  slope  of  the  unflapped  airfoil  (from  Reference  1). 
a 

CEXT  is  extended  chord  and  is  defined  as  Cj  for  the  single -slotted  flap 
and  c  for  the  double-slotted  flap  (from  Figure  4). 

2.1.4  Maximum  Lift  Increment 

The  maximum  lift  increments  obtained  from  thin  airfoil  theory  for  flapped  airfoils 
are  strongly  dependent  on  the  flow  characteristics  of  the  unflapped  airfoil  at  maxi¬ 
mum  lift.  The  plienomena  for  conventional  airfoils  are  directly  related  to  the  air¬ 
foil  leading  edge  radius  and  trailing  edge  angle.  These  parameters  are  determined 
by  the  airfoil  family  and  thickness  distribution  (see  References  1  and  2). 

Basically,  the  deflection  of  the  trailing  edge  introduces  camber  to  the  airfoil  and 
results  in  a  lift  increase.  Thin  airfoil  theory  indicates  a  basic  load  distribution 
from  angle  of  attack  and  an  incremental  load  distribution  due  to  camber  or  flap 
deflection,  as  shown  in  the  following  sketch. 


Basic  Load  Distribution  Load  Distribution  Due  to  Flap  Deflection 
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The  loading  caused  by  flap  deflection  has  a  peak  at  the  flap  hinge  line  that  diminishes 
to  zero  near  the  flap  trailing  edge  and  reduces  to  a  much  lower  v^lue  at  the  leading 
edge.  Nominally,  the  forward  pressure  gradient  is  relieved  and  me  gradient  over 
the  flap  is  increased,  causing  separation  even  at  low  flap  deflection.  This  trend  is 
noticeable  in  Figure  5,  which  indicates  the  turning  efficiency  of  plain  flaps  at  0. 
The  efficiency  is  reduced  by  separation  and  begins  to  fall  off  at  flap  deflections  above 
15  degrees.  The  flow  separation  is  confined  aft  of  the  hinge  line  and  does  not  spread 
forward  until  the  flap  deflection  becomes  increasingly  large  or  the  flapped  airfoil 
reaches  a  maximum  angle  of  attack  and  stalls.  Maximum  lift  can  then  be  defined  as 
the  lift  limited  by  either  flow  separation  that  progresses  forward  of  the  flap  or  flow 
separation  from  the  leading  edge. 

The  theoretical  relationship  from  thin  airfoil  theory  was  applied  to  the  available 
maximum  lift  data  on  trailing  edge  flaps.  The  theoretical  ratio  of  the  maximum 
lift  increment  to  the  lift  increment  developed  at  a  =  0  is  based  on  a  criterion  of 
leading  edge  separation  on  thin  airfoils.  As  a  consequence,  empirical  factors 
were  required  to  correlate  the  experimental  data  on  airfoils  with  finite  thickness. 

The  correlation  parameters  account  for  leading  edge  radius,  thickness,  flap  chord, 
and  flap  deflection. 


The  resulting  expression  for  maximum  lift  increment  for  plain,  single-  and  double - 
slotted  flaps  uses  the  lift  increment  at  #  =0  from  Section  2.1.2  and  is  given  by  the 
following  expression: 


.Ac 

/  I 


Ac 


=  Ac 


nay.  ’ 


max 


.  Ac 

a-°\  Vo/ 


K  K 
T  6 


(27) 


TH 


where  Ac.  is  the  predicted  lift  increment  (from  Section  2. 1.2). 

Xj  a 

a  =  0 


.Ac 


max 


Ac. 


is  the  theoretical  relationship  accounting  for  flap  chord 


«-0  L 


TH 

(from  Figure  2). 


is  the  empirical  factor  developed  from  experimental  data  for  a  flap  at 
an  optimum  deflection  angle  (from  Figures  7  and  8) . 


Kg  is  tha  factor  accounting  for  changes  in  flap  deflection  from  the  optimum 
deflection  (from  Figure  9). 
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Maximu.Ti  Lift  Correlation  Factor 


Figure  8.  Maximum  Lift  Correlation  Factor  for  Double -Slotted  Trailing 
Edge  Flaps  Versus  Leading  Edge  Radius/Thickness  Ratio. 
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The  factor  K™  was  developed  to  account  for  flow  separation  from  either  the  leading 
or  trailing  edge  on  airfoils  with  finite  thickness. 

'The  parameter,  leading  edge radius/thickness ratio,  shown  in  Figures  7and  Sis  based 
on  leading  edge  radius  and  thickness  ratio  of  the  basic  airfoil.  These  terms  corre¬ 
late  the  data  and  give  an  adequate  description  of  maximum  lift  increment  for  the 
trailing  edge  flaps  considered.  The  data  trends  from  recent  General  Dynamics  two- 
dimensional  tests  are  also  shown  in  Figures  7  and  9.  These  tests  were  conducted 
in  a  facility  that  uses  blowing  slots  on  the  walls  to  reduce  separation  effects. 

Prediction  of  the  maximum  lift  of  single -slotted  flaps  with  full  chord  extension  re¬ 
quires  an  additional  factor.  This  turm,  K^,,  is  shown  in  Figure  9. 

2.1.5  Pitching  Moment  Increment 

The  equations  developed  from  thin  airfoil  theory,  Section  2.1.1,  show  that  the  pitch¬ 
ing  moment  increment  due  to  flap  deflection  is 


The  moment  was  then  directly  related  to  the  lift  increment  at  a  =  0  by  the  following 
expression: 


Ac 

m 


c 

l 


6 


The  resulting  values  for  the  theoretical  term  are  given  in  Figure  3.  These  values 
are  used  in  conjunction  with  the  predicted  lift  increment  at  a  =0  to  develop  empiri¬ 
cal  factors  and  correlate  experimental  data.  The  expression  for  the  pitching 
moment  increment  becomes 


Ac 

m 


where  Ac^  is  the  predicted  lift  increment  (from  Section  2. 1.2). 

a=0 


(28) 
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Flap  Angle  Correction  Fact 
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Figure  9.  Flap  Angle  Correction  Factor  Versus  Flap  Deflection. 
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is  the  theoretical  relationship  accounting  for  flap  chord  (from 


Figure  3). 


Kyjj  is  the  empirical  factor  developed  from  experimental  data  (from  Figure 

10). 


Thin  airfoil  theory  gives  a  good  approximation  of  moment  for  trailing  edge  flaps  for 
which  the  flow  remains  attached.  However,  as  flap  chord  and  thickness  are  in¬ 
creased,  factors  are  required  to  account  for  flow  separation  as  shown  in  Figure  10. 

2.1.6  Pitching  Moment  Curve  Slope 

Deflection  of  a  trailing  edge  flap  on  a  thin  airfoil  introduces  camber  to  the  basic 
airfoil.  The  thin  airfoil  theory  discussed  in  Section  2. 1. 1  indicates  that  camber 
does  not  change  the  moment  curve  slope  of  the  airfoil.  Experimental  data  verify 
that  there  is  no  change  in  slope  for  the  angle  of  attack  and  flap  deflection  ranges 
of  plain  flaps  for  which  flow  is  attached.  At  angles  of  attack  above  the  linear  range, 
the  flow  over  the  flap  separates  and  the  additional  loading  from  the  flap  is  lost  and 
pitch-up  results.  Below  the  linear  range,  the  flow  separates  on  the  underside  of  the 
airfoil,  causing  a  forward  shift  in  center  of  pressure  and  a  resulting  nose-up 
moment  change. 


This  investigation  only  concerns  itself  with  the  linear  variation  of  the  pitching  mo¬ 
ment  curve  slope.  For  plain  flaps,  the  assumption  is  that  there  is  ho  change  in 
slope.  Trailing  edge  flaps  with  translating  motion  cause  a  change  in  moment  curve 
slope.  This  change  is  related  directly  to  chord  extension  as  indicated  by  References 
19,  22,  and  27,  and  is  described  by  the  following  expression: 


Ac 

m 


-0.25  K 

ac 


(29) 


where  K_ „  is  an  empirical  factor  to  account  for  chord  extension  (from  Figure  11). 

aC 

CEXT  *s  extended  chord  which  is  defined  as  c^^  for  the  single -slotted 
flap  and  C£  for  the  double -slotted  flap  (from  Figure  4). 


The  linear  range  on  the  moment  curve  is  reduced  at  higher  flap  deflections. 
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Drag  Increment 


2.1.7 


The  approach  utilised  for  predicting  the  drag  of  trailing  edge  flaps  is  mainly  em¬ 
pirical.  All  available  data  utilized  in  Section  2.1,1  for  predicting  the  increment 
of  lift  at  a  =  0  was  used  in  the  correlation.  The  correlation  parameter  selected 
was  the  value  of  Ac^a_0  so  that  any  variation  in  drag  with  flap  deflection  was  re¬ 
moved.  This  assumes  that  the  drag  is  directly  related  to  the  lifting  capability  of 
the  flap  and  indirectly  related  to  flap  deflection  and  flap  type . 

Flap  drag  increments  at  a  =  0  for  plain,  single-  and  double -slotted  flaps  are  shown 
in  Figures  12,  13,  and  14. 

2.2  LEADING  EDGE  DEVICES 


2.2.1  Theory 

Thin  airfoil  theory  for  trailing  edge  flaps  is  covered  in  Section  2.1.1.  A  further 
extension  of  this  theory  tc  cover  leading  edge  devices  is  covered  in  this  section. 
The  aerodynamic  lift  and  moment  of  leading  edge  flaps  can  be  obtained  from  trail¬ 
ing  edge  flap  theory  with  a  geometrical  transformation.  Consider  the  flapped  air¬ 
foil  shown  below: 


Geometry  for  Airfoil  With 
a  Trailing  Edge  Flap 


Geometry  for  Airfoil  With 
a  Leading  Edge  Flap 
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Increments  for  Double-Slotted  Flaps  Versus  Lift  Increment  at  a  =  0,  Contd 


Thne  two  flapped  airfoils  have  identical  profiles  to  the  relative  velocity.  The 
following  relationship  may  then  be  written: 


ct'  *  a  +  6f  (30) 

6/  =  -6f  (31) 


Thin  airfoil  theory  for  a  trailing  edge  flap  (Section  2.1.1)  indicates  that  the  lift 
coefficient  is  defined  by  Equation  (6)  as 


Substituting  the  terms  for  the  leading  edge  flap,  we  have 

=  2 v  [ce*  +  (1  +  0^)] 

Differentiating  with  respect  to  6'  at  constant  section  lift, 


Differentiating  with  respect  to  6^  at  constant  a  > 


It  can  also  be  shown  that  at  constant  section  lift, 


(33) 


(34) 


(35) 


and 


(36) 


(37) 


The  leading  edge  flap  lifting  effectiveness  parameter  is  then  computed  using 
Equations  (32),  (34),  and  (35),  and  is  shown  in  Figure  15. 


33 


Referring  back  to  Figure  5,  it  can  be  seen  that  the  lifting  effectiveness  of  a  load¬ 
ing  edge  flap  is  much  less  than  that  of  a  trailing  edge  flap  of  the  same  chord.  The 
maximum  lift  and  moment  paramete  rs  are  shown  in  Figures  16  and  17,  respectively. 


2.2,2  Lift  Increment  at  q  =  0  Degreea 

Leading  edge  devices  change  the  lift  of  the  unflapped  airfoil  by  changing  the  effec¬ 
tive  angle  of  attack.  Unlike  trailing  edge  flaps,  the  deflection  of  a  nose  flap 
causes  a  loss  in  lift  rather  than  an  increase  at  a  =  0.  The  effect  of  these  devices 
may  be  determined  with  the  following  expression: 


where 


Ac  =  c'  6t  _ 

1  -n  LE 

a=0  6 


(38) 


ci  is  the  theoretical  lifting  effectiveness  for  leading  edge  flaps  (from 
*6 

Figure  15). 

6t  „  is  the  deflection  angle  of  the  leading  edge  device. 


The  above  equation  is  used  for  nose  flaps.  The  equation  is  modified  to  account 
for  chord  extensions  associated  with  Krueger  flaps  and  slats  as  shown  below: 


and 


(39) 


(40) 


where 


is  the  ratio  of  extended  chord  to  airfoil  chord  (from  Figure  18). 


2.2.3  Section  Lift  Curve  Slope  Increment 


The  lift  curve  slope  of  the  leading  edge  flap  is  assumed  to  be  the  same  as  that  of 
the  unflapped  airfoil.  For  Krueger  flaps  or  slats,  the  change  in  lift  curve  slope 
due  to  chord  extension  is  given  by 


AC* 

a 


“airfoil 


(41) 
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where  c «  if)  tho  lift  curve  slope  of  the  unflapped  airfoil  (from  Reference  1). 
a. 

c  la  the  extended  chord  (from  Figure  18). 

EX  T 

2.2.4  Maximum  Lift  Increment 

Leading  edge  devices  Increase  the  maximum  lift  of  airfoils  by  lowering  the  nega¬ 
tive  pressures  (suction)  near  the  nose  and  delaying  leading  edge  flow  separation. 
The  maximum  lift  on  the  airfoil  is  achieved  when  the  pressure  on  the  nose  device 
approximates  that  of  the  airfoil  prior  to  stall.  Leading  edge  devices  delay  flow 
separation  at  the  nose  without  significantly  affecting  trailing  edge  flow  and  are 
most  effective  on  thin  airfoils. 


The  relationship  from  the  theory  developed  in  Section  2.2.1  was  applied  to  the 
available  maximum  lift  data  on  leading  edge  devices.  The  theoretical  value  of 
maximum  lift  increment  for  leading  edge  devices  on  a  thin  airfoil  from  Equation 
(37)  is 


(  'max)th 


'c-  -2*  +  2e\ 

V  6  / 
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Empirical  factors  were  developed  to  correlate  all  available  data  on  leading  edge 
devices  on  airfoils  with  finite  thickness.  The  correlation  parameter  from  Section 
2. 1.4,  leading  edge  radius/thickness  ratio,  was  again  used  to  account  for  nose 
shape  and  thickness.  The  maximum  lift  increment  at  optimum  deflection  was  de¬ 
termined  and  then  related  to  the  theoretical  values  shown  in  Figure  16. 


The  following  expression  resulted: 
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where  is  the  maximum  lifting  efficiency  of  leading  edge  devices  (from 

Figure  19)  (data  from  References  34  through  39). 

rjg  is  a  factor  accounting  for  changes  in  flap  deflection  from  the  optimum 
deflection  (from  Figure  20). 

cj  is  the  theoretical  maximum  lifting  effectiveness  (from  Figure  16). 
6  max 
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Figure  18.  Geometry  for  Leading  Edge  Devices. 
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*7 MAX  ~  Maximum  Lifting  Efficiency  for  Leading  Edge  Devices 


Figure  19.  Maximum  lift  Efficiency  for  Leading  Edge  Devices. 
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Deflection  Angle  Correction  Factor 


WT 


Figure  20.  Deflection  Angle  Correction  Factor. 
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is  the  ratio  of  extended  chord  to  airfoil  chord  (from  Figure  18). 


(¥) 

6,  _  is  the  deflection  angle  of  the  leading  edge  device. 
LE 

2.2.5  Pitching  Moment  Increment 

Equation  (36)  shows  that  the  pitching  moment  for  a  nose  flap  is 


m 


=  -c 


m. 


The  c 


term  is  defined  by  the  nose  flap  chord  which  is  equal  to  (1  -  cf/c).  The 
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m  versus  nose  flap  chord  are  shown  in  Figure  17. 
luo 


theoretical  values  of  c 
The  moment  increment  due  to  nose  flap  deflection  then  becomes 
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where  c '  is  the  theoretical  momenv  effectiveness  for  leading  edge  devices 

III  £ 

,  0 

(from  Figure  17). 


This  expression  is  used  for  all  leading  edge  devices. 

2.2.6  Pitching  Moment  Curve  Slope 

As  discussed  in  Section  2. 1. 6,  the  introduction  of  camber  by  either  nose  or  trail¬ 
ing  edge  flap  deflection  causes  no  change  in  the  moment  curve  slope.  Leading 
edge  devices  that  have  translational  motion  produce  a  change  in  slope.  The  follow¬ 
ing  expression  should  be  used  for  slats  and  Kruegers: 


Ac 

TT  a  0'75(cext‘c)  (44) 

i 

where  c  is  the  extended  chord  as  defined  in  Figure  18. 

EXT 

2.2.7  Drag  Increment 

There  is  insufficient  experimental  data  to  empirically  determine  the  effect  of  lead¬ 
ing  edge  devices  on  drag.  Quantitatively  speaking,  leading  edge  devices  cause  an 
upward  shift  in  minimum  drag  with  small  changes  in  drag  level.  For  the  purposes  of 
this  study,  it  is  assumed  that  leading  edge  devices  cause  no  change  in  drag  level. 
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3.0 


PROFILE  DRAG  OF  FUSELAGE,  NACELLES, 
WING  AND  TAIL  SURFACES 


I 


The  profile  drag  is  predicted  by  a  build-up  method  that  calculates  the  flat  plate 
skin  friction  drag  of  each  component  and  accounts  for  three-dimensional  effects 
(interference  and  roughness).  This  approach,  mainly  empirical,  is  more  accu¬ 
rate  than  wind  tunnel  test  data  because  of  Reynolds  number  effects  and  other  items 
which  normally  cannot  be  represented  on  a  wind  tunnel  model. 


3. 1  Wetted  Area 

The  profile  drag  estimation  is  directly  related  to  the  ability  to  determine  wetted 
areas.  The  wetted  area  of  the  wing  and  tail  surfaces  is  readily  obtained  from  the 
planform  area  and  a  factor  to  account  for  thickness: 


SWET  =  KF  X  <planform  area)  (45) 

where  Kp  is  the  planform  area  factor  (from  Figure  21). 

The  fuselage  and  nacelle  wetted  areas  are  determined  by  simplified  formulas  or 
solved  graphically.  One  approach  is  to  measure  the  top  and  side  planview  area, 
and  to  measure  the  circumference,  height  and  width  at  the  constant  section.  The 
wetted  area  is  given  by 

„  CIRCUMFERENCE  „ 

SWET  ”  H  +  W  (STOP  +  SSIDE  <46) 

B  B 


3.2  Profile  Drag 

The  profile  drag  of  each  component  is  defined  as  the  flat  plate  skin  friction  drag 
plus  three-dimensional,  interference,  and  roughness  effects.  The  skin  friction 
drag  is  equal  to  the  wetted  area  multiplied  by  the  skin  friction  coefficient.  The 
coefficient  with  fully  turbulent  flow  is  based  on  Reynolds  number  and  roughness 
effects  and  is  given  in  Figure  22.  A  roughness  factor  value  of  1.2  is  considered 
representative  for  propeller-driven  V/STOL  aircraft  and  is  used  to  determine 
the  value  of  flat  plate  skin  friction  coefficient  from  Figure  22 . 

The  three-dimensional  profile  drag  for  wing,  horizontal,  and  vertical  tail  surfaces 
is  determined  from  the  following  equation  from  Reference  40: 

^f  X  ^WET  4 

CD  =  - ~~ -  [1  +  Kg(t/c)  +  Kp(t/c)  ]  (47) 

°WING 
OR  TAIL 


e  21.  Thickness  Correction  to  Wetted  Area. 


Figure  22.  Turbulent  Skin  Friction  Coefficient  Versus  Reynolds  Number 


where  K  =  supervelocity  term  from  Figure  23. 

s 

Kp  =  60  for  airfoils  with  max  t/c  located  at  30  percent  chord. 

Kp  =  100  for  airfoils  with  max  t/c  located  at  40  to  50  percent  chord. 

The  profile  drag  of  three-dimensional  streamlined  bodies  of  revolution  is  calcu¬ 
lated  with  the  following  equation  from  Reference  40 : 

CD  =  E1  +  0.00l(i/DB)  +  1.  5(Db/*)3/2  +  7(DB/t)G  ]  (48) 

°BODY 

The  second  term  in  the  above  equation  accounts  for  the  boundary  layer  thickness 
increase  for  wrapping  a  flat  plate  into  a  cylinder.  The  third  term  accounts  for 
supervelocity  and  the  fourth  term  pressure  drag.  The  equation  is  used  for  fuse¬ 
lages  and  nacelles.  Deviations  from  the  streamline  shape  should  be  accounted  for 
in  the  fineness  ratio  term  {i/D B)  determined  from  the  length  and  equivalent 
diameter. 


where  S  is  the  frontal  area. 
B 


(49) 


On  nacelles  or  fuselages  where  the  aft  end  is  truncated,  the  fineness  ratio  is  de¬ 
termined  from  an  effective  length  measured  to  a  faired  point. 


Nacelle  Geometry 
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Location  of  Max  Thickness  -  % 


Figure  23.  3uperveloeity  Factor  for 
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The  drag  of  V/STOL  aircraft  fuselages  with  rear  loading  ra  .nps  is  not  adequately 
predicted  by  Equation  (48).  Reference  41  shows  correlations  of  these  type  fuse¬ 
lages  which  resulted  in  new  terms  for  the  fuselage  drag.  These  terms  replace  the 
pressure  drag  term  7{Dg/jt)3  in  Equation  (48)  which  is  valid  for  streamlined 
bodies.  The  resulting  equations  are  indicated  below: 

^fX^WFT  a/o 

CD  = - g - [l+0.001(A/DB)  +  1.5(DB/jt)  ] 

°FTJSE  LAGE 

^[[6(DB/A)5/2-l]  +  5.2(ya/DB)  [1.4  -  (D^^)4]} 

(50) 

For  fuselages  with  afterbodies  that  contract  laterally  as  well  as  longitudinally,  the 
last  term  in  the  equation  drops  out. 
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Fuselage  Geometry 


3.3  Interference  Drag 


The  interference  drag  between  components  is  estimated  with  the  following  equation 
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where  accounts  for  wing -body  interference  (from  Figure  24). 

K2  accounts  for  horizontal  tail-body  interference  (from  Figure  24). 
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Figure  24.  Wing-Body  or  Tail-Body  Interference  Drag  Factors 


Figure  25.  Nacelle-Wing  Interference  Drag  Factors. 


IC  accounts  for  position  of  nacelle  nose  with  respect  to  leading  edge  of 
wing  (from  Figure  25). 


accounts  for  vertical  position  of  the  nacelle  (from  Figure  25). 

If  the  nacelle  is  on  the  wing  tip,  the  interference  is  taken  to  be  one-half  the 
used  for  inboard  nacelles. 

The  drag  of  t^ree -dimensional  bodies  at  angle  of  attack  is  given  by 
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where  S  is  the  frontal  area  of  the  fuselage. 
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value 


(52) 


4.0  SPAN  LOAD  PROGRAM 


The  procedures  developed  in  Sections  2  and  3  are  used  with  the  basic  airfoil  char¬ 
acteristics  in  a  span  load  digital  computer  program.  This  program  was  written 
at  Convair  and  uses  linear  and  nonlinear  flapped  airfoil  section  data  to  predict  the 
three-dimensional  flapped  wing  characteristics.  The  span  load  program,  original¬ 
ly  set  up  to  use  two-dimensional  experimental  data  from  the  Convair  wind  tunnel, 
has  been  modified  to  use  the  predicted  characteristics  for  this  study. 

The  program  is  based  on  the  modified  Weissinger  L-method  described  by  Gray  and 
Schenk  (Reference  42)  and  by  Holt  (Reference  43).  The  method  represents  the  wing 
with  a  system  of  horseshoe  vortices  extending  aft  from  the  quarter-chord  line. 

The  strengths  of  the  vortices,  and  the  corresponding  load  distribution,  are  ex¬ 
pressed  in  terms  of  the  downwash  angles  at  the  three -quarter -chord  points  for  a 
set  of  spanwise  control  stations.  Using  one  control  station  for  each  horseshoe 
vortex  leads  to  a  set  of  simultaneous  algebraic  equations  which  can  be  solved  for 
the  vortex  strengths  required  to  produce  a  specified  set  of  downwash  angles.  The 
coefficients  in  these  equations  are  calculated  from  the  geometry  of  the  wing  plan- 
form.  The  planform  parameters  required  by  the  program  are  the  wing  span,  as¬ 
pect  ratio,  taper  ratio,  and  quarter -chord -line  sweep  angle. 

This  basic  span  load  distribution  program  is  used  in  an  iterative  procedure  to 
calculate  the  spanwise  variation  of  lift  satisfying  angle  of  attack  relationships  at 
each  wing  section,  as  illustrated  in  the  following  sketch. 


Airfoil  Angles  of  Attack 


Airfoil  Lift  Curve 
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Nonlinear  section  lift  curves  (two-dimensional  data)  are  specified  for  each  span- 
wise  control  station,  and  the  spanwise  variation  of  geometric  angle  of  attack  is 
known  from  basic  wing  geometry.  The  basic  set  of  algebraic  equations  represent¬ 
ing  the  horseshoe  vortex  simulation  of  the  lifting  wing  is  used  to  extract  the  induced 
angles  for  an  assumed  initial  load  distribution  for  equality.  This  process  is  cycled 
through  a  controlled  iterative  procedure  until  the  equality  is  satisfied  within  an 
acceptable  tolerance. 

The  span  load  program  in  essence  accepts  the  two-dimensional  lift  curves  and  pre¬ 
dicts  not  only  the  linear  region  of  the  three-dimensional  lift  curve  but  also  the  stall 
characteristics.  The  program  accepts  as  additional  input  two-dimensional  drag 
and  pitching  moment  data.  From  these  data  and  the  calculated  local  lift  coefficients, 
the  program  determines  the  spanwise  moment  and  profile  drag  distributions.  The 
local  drag  and  pitching  moment  are  integrated  numerically  to  obtain  the  coefficients 
for  segments  of  the  wing  which  are  subsequently  used  to  determine  the  lift,  longi¬ 
tudinal  force  and  moment  coefficients  of  a  wing  immersed  in  a  slipstream. 
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5.0  DEVELOPMENT  OF  LIFT  AND  LONGITUDINAL  FORCES  ON 
A  WING-NACELLE -FLAP  COMBINATION  IMMERSED  IN  A 
PROPELLER  SLIPSTREAM 


Due  to  the  limitations  imposed  by  the  existing  semiempirical  and  theoretical 
methods  for  estimating  the  lift  and  longitudinal  force  coefficients  for  a  wing 
partially  immersed  in  a  propeller  slipstream,  a  program  was  initiated  at  Con- 
vair  to  investigate  a  more  basic  approach  to  the  problem.  This  approach  depends 
on  a  method  of  predicting  two-dimensional  data  and  was,  in  part,  first  conceived 
by  Canadair  Limited  in  the  early  1960's  for  use  on  the  proposed  CL-62  aircraft. 
Although  an  attempt  was  made  in  Reference  44  to  summarize  the  Canadair  ap¬ 
proach,  a  rigorous  analysis  has  never  been  published.  This  section  documents, 
in  detail,  the  analytical  reasoning  behind  the  method  and  describes  the  improve¬ 
ments  which  have  been  incorporated  in  it  during  this  effort. 

There  are  basically  four  terms  which  contribute  to  the  lift  and  longitudinal  force 
coefficients  of  an  interacting  wing/propeller  configuration: 

1.  Forces  due  to  the  inclination  of  the  thrust  vector  to  the  free-stream  flow 
(herein  defined  as  direct  thrust  forces) . 

2.  Forces  due  to  the  free-stream  flow  acting  on  the  portion  of  the  wing  out¬ 
side  the  propeller  slipstream  (herein  defined  as  free-stream  forces). 

3.  Forces  due  to  the  superimposed  slipstream/free-stream  flow  acting  on 
the  portion  of  the  wing  inside  the  propeller  slipstream  (herein  defined  as 
slipstream  forces). 

4.  Forces  due  to  the  engine  nacelle  located  in  the  superimposed  slipstream/ 
free-stream  flow  (herein  defined  as  nacelle  forces). 

In  the  present  analysis,  it  is  assumed  that  these  four  items  can  be  analyzed  in¬ 
dependently,  being  functions  only  of  configuration  geometry,  thrust  vector,  T, 
and  thrust  inclination  angle,  a-p.  The  relative  orientation  of  these  quantities  is 
presented  in  Figure  26. 

5. 1  Direct  Thrust  Forces 


It  is  apparent  from  Figure  26  that  the  inclination  of  the  thrust  vector  to  the  free 
stream  results  in  a  direct  contribution  of  this  vector  to  the  lift  and  longitudinal 
force  terms.  The  contributions  are,  respectively, 
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Figure  26.  Propeller  Slipstream  Notation, 


(53) 


where  S  =  reference  wing  area. 

It  should  be  mentioned  here  that  all  coefficients  will  be  referenced  to  slipstream 
dynamic  pressure,  ci  ,  to  eliminate  the  problem  of  infinite  values  in  the  hover 
condition  where  the  free -stream  dynamic  pressure,  q,  goes  to  zero. 


5.2  Free -Stream  Forces 

Since  the  lift  and  drag  characteristics  on  the  free-stream  segment  of  the  wing  are 
assumed  to  be  independent  of  the  slipstream  flow  (and  hence  thrust  coefficient), 
these  forces  may  be  determined  quite  easily  knowing  the  spanwise  variation  of  the 
power-off  section  lift  and  drag  characteristics  of  the  wing  in  question. 
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One  simply  applies  an  integration  technique  to  these  spanwise  variations  over  the 
free-stream  segment  of  the  wing.  Typical  spanwise  variations  of  power-off  lift 
and  drag  coefficients  (c^  and  c^  from  the  span  load  program  described  in  Section 
4.0)  are  shown  in  Figures  27  and  28  respectively  at  an  arbitrary  angle  of  attack. 
The  quantities  (c^c)  and  (c^c)  are  considered  to  be  weighted  in  the  sense  that  the 
section  lift  and  drag  coefficients  are  multiplied  by  the  local  wing  chord,  c,  at  each 
section.  Also  included  on  these  figures  is  a  representative  model  showing  how  the 
flow  over  the  wing  might  be  divided  between  the  free-stream  and  slipstream  re¬ 
gions.  The  shaded  areas  on  both  figures  represent  the  free-stream  segments. 
This  model  will  be  used  in  the  derivation  of  the  expressions  for  determining  the 
lift  and  drag  coefficients  on  the  two  segments  of  the  wing. 

The  quantity  t)  shown  on  Figures  27  and  28  is  a  nondimensional  spanwise  quantity 
defined  as 


1i I  =V  (M) 

2 

where  y  =  spanwise  distance  measured  from  the  wing  root  section, 
b  =  wing  span. 


During  the  course  of  the  analysis,  7)L  will  be  assumed  constant  and  located  at  the 
inboard  tip  of  the  most  inboard  propeller.  In  effect,  this  assumes  that  the  slip¬ 
stream  tube  remains  tubular  rather  than  contracting  with  increasing  thrust  coeffi¬ 
cient,  as  is  normally  the  case.  The  primary  purpose  of  making  this  assumption 
was  to  reduce  computer  time.  Thus,  future  improvement  in  the  method  might 
make  use  of  the  relationship 


where  d  is  the  contracted  stream -tube  diameter. 


The  model  is  further  simplified  by  assuming  that  the  outboard  boundary  of  the  fully 
developed  slipstream  extends  beyond  the  wing  tip  to  eliminate  the  necessity  of  in¬ 
cluding  additional  terms  in  the  free-stream  solution. 

The  integration  of  the  shaded  area  in  Figure  27  and  the  resultant  calculation  of  the 
lift  coefficient  on  the  free-stream  segment  is  carried  out  within  the  operation  of  the 
span  load  program.  The  result  can  be  expressed  mathematically  as 
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This  operation  is  carried  out  at  several  angles  of  attack  throughout  the  entire  lift 
curve.  Using  the  results  obtained  in  the  linear  portion  of  the  lift  curve  yields  the 
following  solution  for  the  lift  curve  slope  of  the  free -stream  segment: 


Aa 


(59) 


where  a  =  wing  angle  of  attack  relative  to  the  root  chord. 


Finally,  the  lift  coefficient  on  the  free-stream  section  at  any  angle  of  attack  a^g 
can  be  found  from  the  relation 

°l =  (°L  )  <0W  m 

FS  \  a/FS 

where 

aFS  ~  aT  +  *T  ~  aoT 

Li 

iT  =  angle  between  thrust  line  and  local  wing  chord  (see  Figure  26). 

aOT  =  angle  between  local  wing  chord  and  wing  zero  lift  line. 


To  determine  the  drag  coefficient  of  the  free-stream  section,  it  is  necessary  to 
express  the  lift  curve  slope  of  Equation  (60)  in  terms  of  an  effective  geometric 
aspect  ratio.  Basically,  this  means  determining  an  aspect  ratio  which,  when  used 
in  conventional  finite  aspect  ratio  theory,  yields  a  three-dimensional  lift  curve 
slope  equal  to  the  slope  found  in  Equation  (59).  From  finite  aspect  ratio  theory, 
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(61) 


where  a  =  =  wing  three-dimensional  lift  curve  slope. 

W  “w 

AR  =  wing  geometric  aspect  ratio. 

a  =  c.  =  wing  two-dimensional  lift  curve  slope, 
o 
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Note  here  that  the  wing  section  is  assumed  to  be  constant  across  the  span  so  that 
the  section  lift  curve  slope  of  the  wing  remains  constant. 


Note  that  Equation  (61)  assumes  an  elliptical  lift  distribution  over  the  wing  span. 
Modification  of  this  equation  for  use  on  the  wing  section  in  the  free  stream  yields 
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(62) 


where  a„„ 
FS 


lift  curve  slope  on  free-stream  segment  of  the  wing. 


Finally,  solving  Equation  (62)  gives  the  expression  for  the  effective  geometric 
aspect  ratio  of  the  free-stream  section  as 


a 


The  drag  coefficient  of  the  free-stream  segment  can  now  be  described  by  the 
equation 


C  +  C 

DP  °i 
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(64) 


where 


the  profile  drag  coefficient  of  the  free-stream  section, 
the  induced  drag  coefficient  of  the  free-stream  section. 


Following  the  same  approach  used  in  the  lift  expression  of  Equation  (58), 
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As  in  the  lift  case,  the  integration  of  the  shaded  area  of  Figure  28  and  resultant 
calculation  of  Cp  is  carried  out  in  the  span  load  program.  Now,  by  refer- 

encing  the  value  of  Cj^g  to  free-stream  area  Spg  rather  than  total  wing  area, 
the  induced  drag  term  can  be  expressed  as 


60 


2 


IT  AR 


(66) 


The  usefulness  of  the  free-stream  effective  geometric  aspect  ratio  term  (ARpg) 
is  readily  apparent.  This  referencing  of  Cj_  to  free-stream  area  is  necessary 
to  remain  consistent  with  the  effective  geometric  aspect  ratio  terminology. 


The  resultant  value  of  Cj^  must  now  be  referenced  back  to  total  wing  area  to 
remain  consistent  with  the profile  drag  term  determined  by  Equation  (65) , 


(67) 


One  additional  correction  must  be  made  to  the  results  of  Equations  (f  0)  and  (64). 
Since  the  two-dimensional  lift  and  drag  coefficients  are  referenced  to  free-stream 
dynamic  pressure  (q),  and  since  all  coefficients  are  to  be  referenced  to  qg,  Equa¬ 
tions  (60)  and  (64)  must  be  modified  to 
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5.3  Slipstream  Forces 
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The  estimation  of  the  aerodynamic  characteristics  of  a  wing  segment  immersed  in 
a  propeller  slipstream  is  by  far  the  most  difficult  area  of  analysis  associated  with 
a  wing/propeller  configuration.  The  basic  difficulty  lies  in  the  fact  that  the  slip¬ 
stream  and  free-stream  flow  characteristics  over  this  segment  are  superimposed 
in  a  manner  which  varies  throughout  transition  from  hover  (Cp  =  1. 0)  to  high¬ 
speed  flight  (Crp  =0.0). 

AS 
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The  procedures  discussed  in  the  following  paragraphs  make  use  of  the  effective 
geometric  aspect  ratio  terminology  for  each  of  these  two  end  conditions  of  the 
flight  regime,  and  suggest  a  total  equivalent  aspect  ratio  concept  for  bridging 
the  gap  to  include  the  intermediate  thrust  coefficients. 

From  Figure  26,  it  is  seen  that  in  the  general  case  of  intermediate  thrust  coeffi¬ 
cients,  the  superposition  of  the  free -stream  and  slipstream  velocities  results  in 
an  effective  slipstream  flow  which  approaches  the  wing  at  an  angle  ass*  This  is 
the  angle  between  the  resultant  slipstream  velocity  vector  and  the  zero  lift  line  of 
the  wing.  In  equation  form. 


ttnn  ~  0  +  im  “  a  (70) 

SS  *  T  o_  '  ' 

L 

where  0  =  the  angle  through  which  the  slipstream  velocity  vector  is  rotated 
due  to  the  superposition  of  the  free -stream  and  slipstream  flows. 


From  geometry  relations,  0  can  be  expressed  as 

-1 


0  =  sin 


C  sinar 


(71) 


In  the  hover  condition  (CT  =  1.0),  where  the  free-stream  velocity  is  zero.  Equa¬ 
tion  (71)  yields  the  result  0  =  0,  so  that  ctgg  =  iT  -  ao^.  Under  these  conditions, 
the  effective  aerodynamic  lift  on  the  slipstream  segment  can  be  approximated  by  a 
term  which  amounts  to  a  turning  of  the  thrust  through  the  angle  agg  =  iT  -aOL* 
This  term  may  be  thought  of  as  a  supplement  to  the  previously  defined  direct 
thrust  term  and  may  be  expressed  as 


Lr  =  NTsinags  (72) 

Assuming  agg  to  be  a  small  angle  (sin  agg  agg)  and  reducing  to  coefficient  form, 
Equation  (72)  becomes 
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But,  in  the  hover  case,  CTg  =  1.0,  so  that 
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Differentiation  of  Equation  (74)  with  respect  to  ^gg  yields 


a 


H 


(75) 


Note  that  the  slipstream  area  Sgg  is  used  as  a  reference  area  to  again  remain 
consistent  with  the  effective  geometric  aspect  ratio  terminology. 


Finally,  the  effective  geometric  aspect  ratio  of  the  slipstream  section  in  hover 
can  be  derived  from  finite  aspect  ratio  theory: 


a 


Considering  the  other  end  condition  of  high-speed  flight  (C<pg  =  0.0)  where  the 
slipstream  velocity  increment  is  negligible  compared  to  the  free  stream,  Equation 
(71)  yields  the  result  <*  =  a-r*  80  t*13*  ass  =  &T  +  *T  “  a0r  *  This  en<*  condM°n 
thus  reduces  to  a  problem  similar  to  that  of  Section  5.2  w6ere  the  wing  is  subjected 
to  flow  coming  at  it  from  the  free-stream  direction.  The  lift  coefficient  on  this 
slipstream  segment  for  the  high-speed  flight  condition  may  now  be  determined  by 
integrating  the  power-off  span  load  distribution  over  the  slipstream  area. 


Mathematically,  the  integration  of  the  span  load  distribution  and  the  resulting  cal 
culation  of  lift  coefficient  on  the  slipstream  section,  Cr  _,  can- be  expressed  by 
the  relation 
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(77) 


As  in  the  free-stream  segment  of  the  wing,  carrying  out  this  span  load  solution  at 
several  angles  of  attack  gives  results  both  in  the  linear  and  nonlinear  portions  of 
the  lift  curve.  The  linear  results  are  again  useable  in  determining  the  lift  curve 
slope  on  the  slipstream  segment  of  the  wing  as  follows: 


(78) 
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This  lift  curve  slope  of  the  slipstream  segment  in  the  high-speed  flight  condition 
is  converted  to  an  effective  geometric  aspect  ratio  as  follows: 


(79) 


There  now  exist  two  effective  geometric  aspect  ratios  ARH  and  AR^g  which  satisfy 
the  end  conditions  of  hover  and  high-speed  flight  respectively.  In  the  transition 
region  between  these  two  extremes,  a  total  equivalent  aspect  ratio  is  assumed  to 
have  the  form 


“eq  ■  a*h  +  £<“ib-abh>  (80) 

where  V  =  free -stream  velocity. 

V  =  resultant  slipstream  velocity, 
b 

However,  since 
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Equation  (80)  may  be  rewritten  as 

areq  =  "Si  +  /T7Ss  <arhs  -  arh> 


(81) 


This  equivalent  aspect  ratio  may  be  converted  to  an  equivalent  lift  curve  slope  for 
the  slipstream  segment,  using  the  relationship 


and  the  total  lift  coefficient  on  the  slipstream  segment  may  be  calculated  from  the 
equation 
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where 
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5.4  Nacelle  Forces 


The  lift  curve  slope  of  the  nacelle  has  been  estimated  to  be  0.0349  per  degree 
based  on  the  cross-sectional  area  of  the  nacelle.  Since  all  coefficients  are  to 
be  referenced  to  total  wing  area  S,  the  following  relation  is  used: 

„  2 
"dn 

CT  =  0.0349  — —  <89) 

L  4S 

aN 

where  D  =  nacelle  diameter. 

N 

The  lift  on  the  nacelle  can  now  be  found  from  the  relationship 


V  ‘  CL 

N  a 


(90) 


N 

since  the  nacelle  is  located  in  the  slipstream  flow. 


The  drag  on  the  nacelle  is  determined  from  the  methods  of  Section  3.0.  This 
value  must  be  corrected  for  thrust  coefficient  by  the  relation 


(i-cTg) 


(91) 


Finally,  since  the  nacelle  lift  and  drag  terms  are  respectively  perpendicular  and 
parallel  to  the  resultant  slipstream  velocity  vector,  they  must  be  put  into  the 
free-stream  coordinate  system.  Thus, 
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CL  COS  (of  -  0) 
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CL  sin  (aT  -  0) 


-  sin  (aT-0) 
+  CD  cos(aT-0) 
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where  N  =  total  number  of  propellers. 
5. 5  Summary  Equations 


Summarizing  the  results  of  the  previous  sections  yields  the  following  total  equa¬ 
tions  for  the  lift  and  longitudinal  forces  on  a  wing/propeller-nacelle  configuration: 


66 


Direct  Thrust 
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The  two-dimensional  stall  characteristics  as  determined  from  Section  2.0  for  the 
airfoil  section  or  high-lift  devices  being  investigated  are  used  in  the  span  load 
program  to  calculate  the  maximum  lift  coefficients  of  the  free -stream  and  slip¬ 
stream  segments  of  the  wing.  Once  these  values  are  reached,  a  flat  plate  theory 
is  applied  which  calculates  the  lift  and  drag  on  the  wing  at  angles  of  attack  past 
stall.  These  equations  for  the  free -stream  forces  are 


(96) 
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where  CN  is  calculated  at 
FS 
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“fs  “max. 


In  the  same  manner,  the  slipstream  forces  are  found  to  be 


C  —  C  cos  q/ 
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C  =  C  T  sin  (v 
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where  C  is  calculated  at 

SS 

c  =  c 
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SS  MAX _ 


“ss  “max. 


These  summary  equations  and  stall  limits  have  been  programmed  using  standard 
FORTRAN  IV  language.  A  program  description,  listing,  and  required  input  are 
presented  in  Reference  45. 
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6.0  MOMENTS  ON  A  WING-NACELLE-FLAP  COMBINATION  IMMERSED 

IN  A  PROPELLER  SLIPSTREAM 


The  methods  currently  available  for  estimating  the  lift  and  drag  characteristics  of 
propeller-wing-flap  V/STOL  configurations  do  not  predict  the  effects  of  propeller 
slipstream  on  pitching  moments.  Methods  are  presented  herein  for  predicting 
power  effects  on  zero-lift  pitching  moments  and  the  pitching  moment  curve  slope 
of  a  wing  immersed  in  a  slipstream.  These  methods,  developed  at  Convair,  are 
described  in  References  46  and  47  where  the  backup  data  correlations  are  shown. 
Use  is  made  of  the  span  load  program  to  obtain  the  basic  power-off  characteris¬ 
tics  of  the  wing-flap  system  being  analyzed,  and  analytical  expressions  are  pre¬ 
sented  to  correct  these  terms  for  power  effects. 

There  are  three  basic  contributions  to  the  moment  of  a  wing  immersed  in  a  pro¬ 
peller  slipstream: 

1.  Forces  due  to  the  basic  wing. 

2.  Forces  due  to  the  inclination  of  the  thrust  vector  (herein  defined  as  direct 
thrust  forces). 

3.  Forces  due  to  the  lift  augmentation. 

This  analysis  assumes  that  these  forces  can  be  analyzed  separately,  being  func¬ 
tions  only  of  configuration  geometry  and  thrust  coefficient. 

6.1  ZERO  LIFT  PITCHING  MOMENT 

The  individual  contributions  to  the  zero  lift  pitching  moment  of  a  flap  immersed  in 
a  propeller  slipstream  are  shown  in  Figure  29  and  described  below. 

6.1.1  Power-Off  Contribution 
The  power-off  contribution  is 

(*Cm  0)  =  (^TsH^mJ  (100) 

\  /po  \  M  /F 

where  /ACm  )  is  the  zero  lift  pitching  moment  increment  of  the  flap, determined 
'F  from  the  span  load  program. 
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Figure  29.  Power  Effects  on  Zero  Lift  Pitching  Moment. 


6.1.2  Direct  Thrust  Contribution 
The  direct  thrust  contribution  is 


(101) 


where 


is  the  zero  lift  pitching  moment  increment  of  that  portion 
of  the  flap  which  is  immersed  in  the  propeller  slips cream, 
determined  by  the  span  load  program. 


The  spanwise  locations  of  the  flap  ends  are  shown  in  Figure  30  and  defined  below: 


*1 


(102) 


*?o  =  1'° 


(103) 
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Figure  30.  Geometry  for  Flaps  in  Propeller  Slipstream. 

As  in  Section  5.0,  r?  will  be  assumed  constant  and  located  at  the  inboard  end  of 
the  propeller. 

There  is  an  additional  thrust  contribution  from  the  geometrical  relationship  be¬ 
tween  the  thrust  line  and  the  moment  reference  center: 


1 


71 


(104) 


where  c  =  mean  aerodynamic  chord. 

Lp  =  moment  arm  from  quarter  chord  of  the  MAC  to  the  propeller  plane. 

6.1.3  Lift  Augmentation  Contribution 

The  lift  augmentation  contribution  is 


(105) 


6.2  PITCHING  MOMENT  CURVE  SLOPE 

The  calculation  of  the  effects  of  propeller  slipstream  on  the  pitching  moment  curve 
slope  resolves  itself  into  determining  power  effects.  The  approach  used  in  the 
method  is  to  estimate  power-off  lift  curve  slope  and  the  direct  thrust  and  lift  aug¬ 
mentation  contributions.  The  power-on  lift  curve  slope  is  resolved  into  the  indi¬ 
vidual  terms  shown  in  Figure  31  and  described  in  the  following  paragraphs. 

6.2.1  Power-Off  Contribution 


The  power-off  pitching  moment  curve  slope  is  corrected  to  slipstream  dynamic 
pressure  in  the  following  manner-. 

,1( 


where 


is  the  pitching  moment  curve  slope  of  the  basic  wing,  deter¬ 
mined  by  the  span  load  program. 


6.2.2  Direct  Thrust  Contribution 


The  direct  thrust  contribution  due  to  propeller  normal  force  and  effective  thrust 
line  displacement  has  been  determined  from  empirical  correlations.  This  vari¬ 
ation  is  shown  in  Figure  32  and  is  assumed  applicable  for  all  configurations. 
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6.2.3  Lift  Augmentation  Contribution 


Since  the  lift  augmentation  is  a  result  of  the  increased  mass  flow  in  the  wiug 
stream  tube  due  to  the  higher  velocity  in  the  propeller  slipstream,  the  pitching 
moments  are  determined  from  the  geometry  of  that  portion  of  the  wing  immersed 
in  the  slipstream. 


where 


(c-ju-(c"»a»w„ 


The  (C~  \  and  /CT  \  terms  are  respectively  the  power-off  pitching  moment 
\  °7SS  \  <*/SS 

curve  slope  and  the  lift  curve  slope  on  the  slipstream  segments  of  the  wing, 
determined  by  the  span  load  program. 
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6.2.4  Nacelle  and  Fuselage  Contributions 


Since  the  nacelles  are  in  the  presence  of  the  slipsti  am,  both  the  power-off 
and  lift  augmentation  contributions  to  the  lift  curve  slope  are  used  to  obtain 
the  nacelle  pitching  moment  curve  slope.  Also,  an  additional  multiplying  term 
has  been  included  because  it  is  expected  that  the  destabilizing  effect  of  the  na¬ 
celles  would  increase  in  direct  proportion  to  the  increase  in  lift  curve  slope 
arising  from  the  lift  augmentation  term.  The  reason  for  including  this  term  is 
intuitive,  but  it  is  needed  to  achieve  correlation,  and  physically,  the  nacelle 
stability  contribution  is  related  to  the  wing  upwash  field  in  the  vicinity  of  the 
nacelle,  which  in  tum  is  related  to  the  aerodynamic  contributions  to  lift  curve 
slope.  The  contribution  of  the  nacelles  is  then  given  by 
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and  N  =  number  of  propellers. 

Kn  =  nacelle  stability  coefficient  (from  Figure  5-16  in  Reference  48). 

Wfj  =  nacelle  width. 

A,  -  nacelle  length. 

N 

c  =  wing  chord  at  nacelle  location. 
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(RAD-1)  for  small  a's. 


=  total  power-on  lift  carve  slope. 


The  pitching  moment  contribution  due  to  the  body  is  found  from  the  relation 


(112) 


kbwb2lb 

Sc 


where  Kg  =  fuselage  stability  coefficient  (from  Figure  5-16  in  Reference  48). 
Wg  =  fuselage  width. 

Lg  =  fuselage  Jengln. 

6.3  SUMMARY  EQUATIONS 


The  power  effects  on  the  total  zero  lift  pitching  moment  coefficient  increments  due 
to  flap  deflection  are  given  by  the  following  equation : 


POWER  OFF  DIRECT  THRUST 


LIFT  AUGMENTATION 

+^^(AHS 


(113) 


These  increments  were  patterned  after  the  relationships  of  Reference  49. 


The  total  pitching  moment  curve  slope  at  any  thrust  coefficient  is  calculated  with 
the  following  equation: 
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The  combination  of  Equations  (113)  and  (114)  provides  values  of  pitching  moment 
for  any  value  of  or  These  equations  and  the  lift  and  longitudinal 

force  equations  have  been  incorporated  into  the  computer  program  described  in 
Reference  45  to  supply  tail-off  aerodynamic  data  to  the  performance  program. 
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7.0  DOWNWASH  CHARACTERISTICS 


The  basic  approach  of  this  procedure  is  to  develop  equations  for  calculating  down- 
wash  velocities  in  the  plane  of  symmetry,  at  the  longitudinal  and  vertical  stations 
of  the  horizontal  tail,  as  given  by  the  Biot-Savart  equation.  The  span  load  distri¬ 
butions  necessary  to  carry  out  the  computations  can  be  obtained  from  the  existing 
computer  program  (described  in  Section  4.0)  based  on  the  method  of  Reference  42. 
Finally,  the  downwash  thus  obtained  is  modified  to  account  for  the  effects  of  a 
slipstream  by  a  procedure  similar  to  that  given  in  Reference  50. 

As  a  starting  point  in  the  analysis,  equations  are  developed  for  calculating  down- 
wash  angle  in  the  plane  of  symmetry  for  any  longitudinal  or  vertical  position  of 
the  horizontal  tail.  An  undeflected,  undistorted  wake  is  assumed.  Since  span 
load  distributions  are  obtained  from  Reference  42,  the  same  conceptual  model  of 
horseshoe  vortices  which  replace  the  circulation  distribution  of  the  actual  wing 
will  be  considered  in  calculating  the  induced  velocities  at  the  tail  plane.  Also, 
the  same  general  form  of  the  Biot-Savart  equation  is  used  for  calculating  the 
downwash  velocity  at  a  point  P  induced  by  an  elemental  vortex  filament: 
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Ar  (cos a -cos 0) 
4^R 
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The  angles  and  distances  in  the  above  equation  are  defined  in  Figure  33.  The  com¬ 
ponent  qp  is  the  total  tangential  velocity  component  induced  at  point  P  by  the  vor¬ 
tex  element,  and  wp  is  the  downwash  velocity  component. 

The  planform  geometry  and  typical  arrangements  of  the  horseshoe  vortices  (cor¬ 
responding  to  those  in  Reference  42)  are  shown  in  Figure  34.  The  tail  length  is 
defined  in  the  usual  sense,  i.e.,  distance  from  the  quarter  chord  of  the  wing 
mean  aerodynamic  chord  to  the  quarter  chord  of  the  horizontal  tail  mean  aero¬ 
dynamic  chord.  In  deriving  the  equations,  the  origin  of  the  x  axis  is  at  the 
quarter  chord  of  the  wing  root  chord,  so  that  the  tail  length  is  defined  by  the 
following  equation: 


4h  =  vKirrWo/i  <U5> 

On  Figure  34,  77.  refers  to  the  spanwise  location  of  the  midpoint  of  the  bound  vor¬ 
tex  for  any  chosen  number  of  spanwise  stations,  and  A77  is  the  bound  vortex  span 
of  each  elemental  horseshoe  vortex  system. 
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Figure  33.  Vortex  Element  Geometry. 


From  the  geometry  for  the  tilt-wing  in  Figure  35,  it  is  seen  that  the  problem  be¬ 
comes  that  of  calculating  the  downwash  at  a  point  aft  of  the  wing  root  quarter  chord 
and  above  the  wing  wake.  These  distances  vary  as  a  function  of  the  wing  tilt  angle, 
angle  of  attack,  and  wake  location.  The  longitudinal  distance  can  be  written  di¬ 
rectly  in  terms  of  tilt  angle,  angle  of  attack,  and  airplane  geometry  as  follows: 

y =[v  KHr)UnAo/4  ■ (I  • cos  ‘w>*piv]  ™(«w  -  ‘w> 

+  (zH  -  ZPIV  -  Xpivsin  iw)sin  (aw  -  iw) 

(116) 


The  vertical  distance  is  given  by 


y=y+ho+h 


(117) 


The  first  two  te  ms  in  the  above  equation  are  dependent  only  on  the  tilt  angle,  angle 
of  attack,  and  airplane  geometry,  as  follows: 
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Figure  34.  P] 


Figure  35.  Side  View  of  Wake  Geometry  for  Tilt-Wing  Configuration. 


ZH;  =  (Zh  "  ZPIV  "  xPIVsin  iw^Cos^aW  “ 
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hQ  =  .75  cr  sin  iw  (119) 

The  last  term  is  the  wake  displacement,  which  requires  an  auxiliary  calculation. 
The  equations  for  calculating  wake  displacement  are  developed  in  subsequent 
paragraphs. 


The  next  step  in  the  analysis  is  to  apply  the  Biot-Savart  equation  to  an  element  of 
the  horseshoe  vortex  system.  With  all  distances  measured  on  semispans,  the  total 
downwash  contribution  of  two  elemental  bound  vortices  at  equivalent  spanwise  sta¬ 
tions  on  the  left  and  right  wing  panels  is  given  by  the  following  equation: 
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Contribution  of  the  trailing  vortices  of  the  left  and  right  wing  panels  is 

(%"-’>ita,lAc/4 

Vfo +f f + (v  -  VanAc/4)2+(V)2  J 

(121) 

Then  the  total  downwash  is  obtained  as  the  summation  of  Equations  (120)  and  (121) 
for  any  specified  number,  i,  of  elemental  horseshoe  vortices.  All  of  the  terms 
in  Equations  (120)and  (121)  are  defined  except  the  distance  from  the  wake  dis¬ 
placement  at  the  longitudinal  station  of  the  horizontal  tail.  The  effects  of  vertical 
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displacement  of  the  wake  are  approximated  by  assuming  that  its  inclination  at  every 
point  is  the  downwash  angle  of  the  wake.  Starting  with  zero  inclination  at  the  wing 
trailing  edge,  the  wake  displacement  is 


[v-v^w-v] 

h  =  J  tan  eodx  (122) 

[xqCOs^w)] 


To  perform  the  wake  displacement  calculations,  Equations  (120)  and  (121)  reduce 
to  the  following  form: 


(123) 

where  x  is  the  distance  behind  the  wing  root  chord  trailing  edge. 

These  equations  determine  the  downwash  in  the  plane  of  symmetry  at  the  longi¬ 
tudinal  and  vertical  location  of  the  quarter  chord  of  the  mean  aerodynamic  chord 
of  the  horizontal  tail.  A  partial  assessment  of  the  validity  of  these  equations  was 
obtained  by  a  comparison  with  some  test  data  from  Reference  51.  The  downwash 
measurements  in  Reference  51  were  obtained  from  vanes  which  were  at  fixed 
locations  in  the  tunnel,  and  therefore  did  not  completely  simulate  a  horizontal 
tail  on  an  actual  configuration.  However,  the  data  do  provide  a  partial  check  and, 
of  more  importance,  provide  information  on  the  effect  of  the  propeller  slipstream 
on  downwash  characteristics.  Since  the  test  data  did  not  represent  a  horizontal 
tail  at  angle  of  attack,  the  comparison  with  calculated  values  was  made  at  zero 
angle  of  attack.  The  comparison  is  shown  m  Figure  36  where  the  downwash  angle 
at  zero  angle  of  attack  and  thrust  coefficient  is  shown  as  a  function  of  vertical 
distance  above  and  below  the  wing  quarter  chord.  As  background  information, 
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the  span  loading  from  which  the  calculated  downwash  angles  were  obtained  is  also 
included  on  Figure  36.  There  is  good  trend  agreement  between  calculated  and 
test  downwash  angles  for  locations  below  the  center  line.  The  test  data  are  notably 
unsymmetrical  about  the  center  line  and  result  in  calculated  downwash  values  that 
are  significantly  higher  than  the  test  values  for  locations  well  above  the  center  line. 
The  agreement  is  considered  to  be  generally  satisfactory  for  trends,  but  adjust¬ 
ments  appear  to  be  required  for  absolute  levels. 

As  noted  in  the  introductory  paragraphs,  power  effects  are  accounted  for  by  modi¬ 
fication  of  the  downwash  calculated  at  zero  thrust  coefficient  by  a  procedure  similar 
to  that  given  in  Reference  42.  In  essence,  the  method  states  that  the  tangent  of  the 
downwash  angle  at  any  thrust  coefficient  is  proportional  to  the  lift  coefficient  at 
that  thrust  coefficient.  The  calculated  downwash  angles  will  thus  have  a  limiting 
value  of  90  degrees  at  a  thrust  coefficient  of  one.  In  equation  form,  the  tangent  of 
the  downwash  angle  at  any  thrust  coefficient  is  simply 

=(tanf)  (CL)  (124) 

S  A  s  s 

Or,  for  the  test  condition  of  Figure  36, 

tan  e  =  .  0449 CL  (125) 

A  comparison  of  the  downwash  calculated  from  the  above  equation  with  the  test  data 
from  Reference  51,  at  zero  angle  of  attack,  is  shown  on  Figure  37.  The  proper 
trend  is  predicted,  but  the  calculated  values  are  higher  than  test  data  throughout 
the  thrust  coefficient  range.  It  was  found  that  improved  correlation  could  be  ob¬ 
tained  by  determining  a  constant  of  proportionality  based  on  an  average  value  ob¬ 
tained  from  all  of  the  available  test  points.  The  equation  becomes 

tan  e  =  .  031 9 (126) 

The  above  equation  is  also  plotted  on  Figure  37.  Constants  of  proportionality  were 
determined  for  test  conditions  in  Reference  51  other  than  those  given  on  Figures 
36  and  37,  and  it  was  found  that  the  average  values  determined  for  all  the  other 
configurations  did  not  differ  greatly  from  .0319.  It  would  therefore  seem  likely 
that  a  constant  of  proportionality  exists  which  would  apply  satisfactorily  to  all 
configurations . 
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Figure  37.  Effect  of  Thrust  Coefficient  on  Downwash  Angle. 
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8.0  WIND  TUNNEL  WALL  CORRECTIONS 


A  study  was  performed  to  investigate  the  validity  and  accuracy  of  wind  tunnel  wall 
corrections  as  applied  to  models  of  tilt-wing  and  deflected-slipstream  V/STOL 
configurations.  It  was  hoped  that  a  significant  amount  of  wind  tunnel  and  flight  test 
data  on  the  same  configuration  would  be  available  to  insure  a  reasonable  conclu¬ 
sion.  However,  a  rigorous  literature  search  revealed  that  useful  flight  test  data 
were  essentially  nonexistent  on  these  types  of  configurations.  In  spite  of  this,  the 
study  was  continued  to  gain  a  basic  familiarity  with  the  available  correction  proce¬ 
dures  in  terms  of  their  relative  complexity  and  the  relative  magnitude  of  the  re¬ 
sulting  corrections. 

A  continuation  of  the  literature  search  indicated  two  basic  correction  procedures 
which  were  considered  worthy  of  further  investigation: 

1.  The  standard  wind  tunnel  wall  correction  procedure  first  proposed  by 
H.  Glauert  in  the  1920's. 

2.  A  more  recent  procedure  developed  by  Mr.  H.  Heyson  of  NASA's 
Langley  Research  Center. 

The  discussion  which  follows  deals  with  the  theoretical  justification  of  each  of  these 
procedures  as  well  as  provides  a  description  of  the  steps  required  in  their  numeri¬ 
cal  calculation.  The  procedures  are  then  applied  to  uncorrected  wind  tunnel  data 
from  models  of  tilt-wing  and  deflected-slipstream  configurations. 

8.1  STANDARD  WALL  CORRECTIONS 

A  complete  description  of  the  standard  wall  correction  theory  is  presented  in 
Reference  52  and  will  not  be  duplicated  herein.  However,  enough  information  will 
be  presented  to  provide  a  general  understanding  of  the  theory. 

8.1.1  Theoretical  Justification 

It  is  well  known  in  fluid  theory  that  a  solid  boundary  near  a  vortex  may  be  simu¬ 
lated  by  the  introduction  of  a  second  vortex  of  opposite  sign  at  some  finite  distance 
away.  This  is  the  basis  for  the  so-called  Method  of  Images  in  which  the  solid 
boundary  corresponds  to  the  zero  streamline  produced  between  the  two  interacting 
vortices.  Since  the  flow  characteristics  about  a  wing  may  be  closely  represented 
mathematically  by  a  simple  vortex  system,  it  is  reasonable  to  assume  that  the 
simulation  of  wind  tunnel  walls  around  a  model  configuration  may  be  appropri¬ 
ately  carried  out  using  this  Method  of  Images. 
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This  is  precisely  the  approach  taken  in  the  standard  wall  correction  procedure. 
Once  the  image  system  is  established,  its  effect  on  the  model  is  theoretically 
identical  to  that  of  the  solid  boundary  it  represents.  Thus,  the  problem  is  re¬ 
duced  to  one  of  finding  the  effect  of  the  image  system  on  the  flow  characteristics 
in  the  neighborhood  of  the  model. 

For  an  image  vortex  located  at  a  distance  from  the  model,  the  induced  upwash 
velocity  at  the  model  as  predicted  by  vortex  theory  becomes 


w  = 


r 

4rrr 


(127) 


In  addition,  the  circulation  for  a  uniformly  loaded  wing  may  be  expressed  as 

HfrW  <i28> 


where  S  =  wing  area 

V  =  f  ree-stream  velocity 
b  =  wing  span 
CL  =  lift  coefficient 


Combining  Equations  (127)  and  (128)  and  expressing  the  distance  in  general  terms 
as  a  constant  times  the  tunnel  height  yields 


w 


\8  tr  kh-jb/ 


(129) 


Now,  rearranging  terms  and  introducing  the  tunnel  width  (B)  such  that  the  product 
of  height  times  width  is  the  tunnel  test  section  area  (C)  gives 


w 

V 


Finally,  defining  the  correction  factor 


6  = 


B 


8  ft  kb 


(130) 


(131) 


and  applying  the  small  angle  assumption  to  Equation  (130),  the  boundary  induced 
upwash  angle  at  the  model  becomes 
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(132) 


w  S  „ 

Aft.  =  tanAai  =  ^=  6^  CL 

Also,  since  the  induced  drag  coefficient  may  be  approximated  by 

CDi  =  CL  (133) 

the  change  in  the  induced  drag  caused  by  the  boundary  induced  upwash  becomes 

ACd.  =  CLAtti  -  6^Cl2  (134) 

Equations  (132)  and  (134)  are  assumed  to  be  applicable  for  all  test  conditions,  even 
though  they  were  strictly  defined  for  the  case  of  a  small  wing  and  uniform  load 
distribution.  This  assumption  is  validated  by  the  fact  that  6  is  varied  to  account 
for  the  type  of  span  load  distribution,  the  ratio  of  model  span  to  tunnel  width,  the 
shape  of  the  test  section,  and  even  the  position  of  the  wing  relative  to  the  tunnel 
center  line.  The  determination  of  5  is  facilitated  by  the  many  charts  presented 
in  Reference  52. 

8.1.2  Calculation  Procedure 

The  use  of  Equations  (132)  and  (134)  for  all  test  conditions  is  made  possible  by 
appropriately  varying  the  correction  factor  5.  Once  this  factor  is  determined  for 
a  given  set  of  test  conditions,  it  remains  constant,  and  the  calculation  of  the 
boundary  induced  upwash  angles  and  resulting  changes  in  induced  drag  at  various 
values  of  lift  coefficient  becomes  a  straightforward  procedure.  Therefore,  be¬ 
cause  of  the  importance  of  this  5  term,  it  is  instructive  to  understand  its  deri¬ 
vation  using  the  charts  presented  in  Reference  52. 

A  quick  examination  of  Reference  52  indicates  the  existence  of  charts  covering 
only  elliptically  or  uniformly  loaded  wings.  Since  most  wings  exhibit  neither  of 
these  loading  conditions,  it  would  appear,  at  first,  that  the  whole  theory  is  rather 
academic.  Closer  examh.ation,  however,  reveals  that  by  using  an  effective  vor¬ 
tex  span  rather  than  the  tctal  wing  span,  the  use  of  the  uniform  load  distribution 
curves  is  quite  proper.  This  is  because  the  shed  vortices  of  an  arbitrary  planform 
wing  rapidly  roll  up  into  a  single  pair  of  vortices  which  exactly  duplicate  the  trail¬ 
ing  vortex  pattern  of  a  uniformly  loaded  wing.  Thus,  most  planforms  to  be  studied 
can  take  advantage  of  the  uniform  loading  charts  of  Reference  52  as  long  as  the 
effective  vortex  span  notation  is  used. 

Given  t’  e  taper  ratio  and  aspect  ratio  of  any  wing  in  question,  Figure  G;25  of  Ref¬ 
erence  52  yields  the  value  of  by/b,  where  by  is  the  vortex  span  downstream  of  the 
model.  The  effective  vortex  span  is  then  determined  from  the  relation 
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(135) 


b  +  b 


Now,  calculating  the  ratio  of  tunnel  height  to  tunnel  width 

hrp 

»T-B-  <136’ 

along  with  the  parameter 

be 

k=—  (137) 

and  assuming  the  model  to  be  located  in  the  center  of  a  closed  rectangular  tunnel, 
the  value  of  6  can  be  found  directly  from  Figure  6:32  of  Reference  52. 

A  similar  procedure  may  be  followed  for  the  case  of  an  open  rectangular  tunnel, 
an  open  and  closed  circular  tunnel,  an  open  and  closed  square  tunnel,  an  open  and 
closed  elliptic  tunnel,  and  an  open  and  closed  circular  arc  tunnel.  In  addition, 
corrections  are  also  presented  for  conditions  where  the  model  is  not  located  in 
the  ce  v’"  1  these  various  tunnels. 

Basically,  the  application  of  the  standard  wall  correction  theory  to  wind  tunnel 
data  is  a  straightforward  task,  the  greatest  problem  being  the  determination  of  the 
factor  6.  Results  of  sample  calculations  on  various  sets  of  wind  tunnel  data  are 
presented  in  Section  8.3. 

8.2  HEYSON'S  CORRECTIONS 

A  complete  description  of  the  basic  Heyson  wall  correction  theory  is  presented  in 
Reference  53.  Only  that  information  deemed  pertinent  to  a  general  understanding 
of  the  theory  is  presented  here. 

8.2.1  Theoretical  Justification 

The  theoretical  justification  of  the  Heyson  method  depends  on  the  application  of  the 
Method  of  Images  to  the  mathematical  representation  of  the  model  wake  by  a  uni¬ 
form  distribution  of  point  doublets  located  along  the  wake  axis.  Note  that  in  th" 
standard  procedure  of  Section  8. 1.1,  the  wake  was  represented  mathematically 
by  a  vortex  system  which  was  assumed  to  trail  horizontally  down  the  tunnel, 
neglecting  any  interference  with  the  tunnel  floor.  Heyson 's  method  accounts  for 
the  downward  deflection  of  the  model  wake  through  the  superposition  of  both  ver¬ 
tical  and  longitudinal  doublets*  The  magnitude  of  this  deflection  is  defined  by  a 


89 


so-called  wake  skew  angle  which,  as  shown  in  Figure  38,  is  the  compliment 
of  the  more  commonly  used  wake  deflection  angle  0 . 


Since  both  longitudinal  and  vertical  velocities  are  induced  from  each  of  these  types 
of  doublets,  four  separate  interference  factors  result  for  each  wake  skew  angle: 
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=  the  interference  factor  associated  with  the  vertical 
velocity  induced  by  the  vertical  doublet 


=  the  interference  factor  associated  with  the  longi¬ 
tudinal  velocity  induced  by  the  vertical  doublet 


=  the  interference  factor  associated  with  the  vertical 
velocity  induced  by  the  longitudinal  doublet 


=  the  interference  factor  associated  wite  the  ongi- 
tudinal  velocity  induced  by  the  longitudinal  doublet 


where 


=  wind  tunnel  cross-sectional  area  (same 

=  =  model  momentum  area. 

4 


as  C  in  the  standard  method). 


The  terms  wQ  and  uQ  are  the  vertical  and  longitudinal  induced  velocities  respec¬ 
tively,  as  determined  from  momentum  theory. 

Heyson's  original  theory,  as  defined  in  Reference  53,  assumed  the  model  to  be 
vanishingly  small  so  that  only  one  line  of  doublets  was  necessary  to  accurately 
represent  the  model  wake.  He  has  since  extended  the  theory  to  cover  several 
specific  model  configurations  including  the  finite  span  swept  wing.  The  method  of 
determining  the  interference  factors  for  these  various  configurations  is  discussed 
in  detail  in  Reference  54.  For  the  case  of  the  finite  span  swept  wing,  the  basic 
procedure  is  to  divide  the  wing  into  several  segments  and  then  represent  the  wake 
of  each  segment  by  its  own  doublet  distribution.  The  resulting  interference  fac¬ 
tors  for  each  segment  are  then  averaged  over  the  model  span  to  determine  a  total 
interference  factor  for  the  wing  in  question. 


Mr.  Heyson  has  made  one  additional  modification  to  his  original  theory  which 
warrants  mention:  the  introduction  of  an  effective  wake  skew  angle  in  place  of 
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Figure  38.  Wake  Skew  Angle  Definition. 


the  old  wake  skew  angle  shown  in  Figure  38.  This  effective  wake  skew  angle  was 
originally  defined  in  Reference  55  as 


v  =  XI?® 

*eff  2 


(138) 


and  was  incorporated  primarily  to  account  for  the  roll-up  of  the  model  wake.  This 
roll-up  results  in  an  effective  wake  d-jflection  angle  which  is  approximately  one-half 
of  that  predicted  at  the  center  of  lift. 

Throughout  the  remainder  of  this  section,  all  interference  factors  calculated  by 
Heyson'.  method  will  be  averaged  values  over  the  model  span  and  will  be  based  on 
the  effective  wake  skew  angle  definition.  A  discussion  of  the  use  of  these  inter¬ 
ference  factors  in  determining  the  final  corrected  lift  and  drag  coefficients  of  an 
arbitrary  model  is  found  in  the  following  section. 
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8.2.2  Calculation  Procedure 


The  general  procedure  for  applying  Heyson's  wind  tunnel  wall  correction  theory  to 
uncorrected  wind  tunnel  data  is  presented  in  Appendix  C  of  Reference  53.  The 
remainder  of  this  section  deals  with  the  basic  equations  contained  in  the  calcula¬ 
tion  procedure  and  includes  a  description  of  the  mathematical  techniques  required 
in  its  computerization. 


The  basic  equations  in  Heyson's  method,  as  given  in  Reference  53,  are  shown  below: 


where  w^ 
L 
D: 


n 


=  reference  velocity 
=  lift  foice 
-  induced  drag  force 
=  tunnel  density 

=  ratio  of  final  induced  velocities  in  the  far  wake  to  the  initial  induced 
velocities  u0  and  wQ  at  the  model  (n  =  2  for  a  wing/propeller  config¬ 
uration). 


A  close  examination  of  Equations  (139),  (140),  and  (141)  revealed  the  possibility  of 
expressing  the  equations  in  a  more  usable  form.  For  example,  with  n  =  2  and  A^  = 
7Tb2 /4 ,  Equation  (139)  may  be  written  as  follows: 
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or 


Now,  writing  Equation  (140)  as 


(143) 


V 


wo  wo 


wh 


(144) 


and  substituting  into  Equation  (141)  yields 


Finally,  for  simplicity,  it  is  normally  assumed  that 


D. 

l 

~L 


D 

L 


so  that  Equation  (145)  may  be  written  as 


(145) 


(146) 


(147) 


Thi  problem  is  thus  reduced  to  one  of  finding  the  values  of  Wq/w^  which  satisfy 
Equation  (147)  for  values  of  V/w^  and  D/L  as  determined  from  Equations  (143) 
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and  (146),  respectively.  Actually,  according  to  Reference  56,  only  one  root  of 
Equation  (147)  is  needed:  that  between  0  and  1.  For  this  case,  it  becomes  possible 
to  use  the  Newton-Raphson  iteration  scheme.  V/wQ  is  then  calculated  from  Equation 
(144)  and  the  value  of  the  wake  skew  angle  is  found  from  the  expression 


X  = 


IT 

180 


arc  tan 


(148) 


which  eliminates  the  problem  of  determining  the  proper  sign  of  the  wake  skew  angle 
Finally,  the  effective  wake  skew  angle  is  determined  from  Equation  (138). 

With  the  proper  value  of  the  effective  skew  angle  known,  it  becomes  necessary  to 
determine  the  four  interference  factors  6w  6w  6U(l»  and  6u  D#  The  pro¬ 
cedure  followed  in  the  computer  program  listed  in  Reference  53  was  to  parametri¬ 
cally  determine  the  interference  factors  at  10-degree  increments  of  xeff  between 
20  and  90  degrees.  The  resulting  values  are  then  used  in  a  slope  iuttrcept  routine 
to  determine  the  correct  values  which  correspond  with  the  value  of  Xeff  • 

The  remainder  of  the  calculation  procedure  is  a  straightforward  task,  and  the  fol¬ 
lowing  equations  resulted: 


A a  =  arc  tan 


(149) 


II 

of  I® 

K)*fsf 

(150) 

Lc  = 

L  cosAa  -  D  sinAa 

(151) 

Dc  = 

L  sinAa  -  D  cosAa 

(152) 

qc 

qc  =  T  q 

(153) 

a  =  a  +Aa 
c 

(154) 

L 

Lc  q  S 

(155) 
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(156) 


;c  q  S 


CT  =  C„  — 


(157) 


A  visual  picture  of  the  effect  of  these  corrections  is  shown  in  Figure  39.  Since  the 

corrected  thrust  coefficient  varies  throughout  a  given  wind  tunnel  run,  as  does 

Cl  and  Cq  ,  some  cross-plotting  must  be  done  to  obtain  a  valid  representation  of 
c  c 

the  corrections  at  a  constant  corrected  thrust  coefficient. 


Resultant 

Force 


Figure  39.  Geometry  for  Heyson's  Wall  Corrections. 
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RESULTS 


8.3 

The  first  data  to  which  the  standard  and  Heyson  wall  corrections  were  applied  are 
presented  in  Reference  56.  In  that  report  a  study  is  described  concerning  an  in¬ 
vestigation  on  the  effects  of  Heyson  *s  wall  corrections  as  applied  to  a  tilt-wing 
model  (without  fuselage)  tested  in  three  different  size  tunnels:  (1)  the  Langley  7- 
x  10-foot  tunnel,  (2)  the  17-foot  diffuser  section  of  this  7-  x  10-foot  tunnel,  and 
(3)  the  Langley  full-scale  30-  x  60-foot  tunnel. 

The  model  tested  had  a  4.23-foot  wing  span  with  two  symmetrically  located  non¬ 
overlapping  propellers,  each  having  a  diameter  of  2  feet.  In  addition,  the  model 
was  equipped  with  a  30-percent-chord  trailing  edge  flap  which  could  be  deflected 
to  an  angle  of  40  degrees.  All  tests  were  run  at  the  same  slipstream  Reynolds 
number.  The  same  sting  attachment  was  used  on  the  model  to  eliminate  these  items 
from  the  correlation. 

The  data  used  in  this  study  were  for  the  flaps -down  configuration,  and  were  obtained 
in  the  7-  x  10-foot  and  30-  x  C9-foot  tunnels  only.  None  of  the  17-foot  diffuser 
section  data  is  presented  because  of  the  questionable  application  of  Heyson 's  cor¬ 
rections  to  these  data.  This  diffuser  section  converges  quite  rapidly  downstream  of 
the  model  and  thus  negates  that  aspect  of  Heyson 's  theory  which  requires  the  doublet 
distribution  to  flow  in  a  straight  line  along  the  floor  following  its  downward  deflec¬ 
tion.  In  addition,  a  close  inspection  of  the  basic  uncorrected  data  showed  the  17-foot 
diffuser  section  data  to  be  inconsistent  with  those  of  the  other  sections. 

Most  of  the  corrections  applied  to  the  data  of  Reference  56  assumed  the  model  to  be 
vanishingly  small,  and  did  not  include  the  new  definition  of  effective  wake  skew  angle. 
This  analysis  corrects  for  both  of  these  deficiencies  in  the  application  of  Heyson 's 
method  to  the  7-  x  10-foot  tunnel  data.  The  30-  X  60-foot  tunnel  data  are  assumed 
to  be  free  of  wall  effects  due  to  the  small  ratio  of  model  to  tunnel  size. 

The  increments  in  induced  upwash  angle  and  induced  drag  coefficient,  as  calculated 
by  the  standard  procedure  of  Section  8. 1.2,  are  applied  to  the  test  data  of  Reference 
56.  The  final  corrected  lift  and  drag  curves  obtained  from  these  standard  correc¬ 
tions  are  shown  in  Figures  40  and  41  respectively. 

The  calculation  of  Heyson ’s  corrections  for  the  data  of  Reference  57  proved  to  be 
considerably  more  difficult  th  .  the  standard  correction  method.  The  basic  cor¬ 
rections  themselves  were  determined  fairly  easily  using  the  computer  program 
listed  in  Reference  54.  However,  a  problem  arose  in  reducing  the  corrected  values 
of  lift,  drag  and  angle  of  attack  to  a  constant  thrust  coefficient  basis.  This  pro¬ 
cedure  was  necessitated  by  the  fact  that  Heyson 's  procedure  corrects  the  tunnel 
dynamic  pressure  and  hence  the  thrust  coefficient  of  every  data  point.  Cross  plots 
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of  the  corrected  values  of  lift,  drag  and  angle  of  attack  were  required.  Data  at 
constant  values  of  corrected  thrust  coefficient  (C^  =4,  8  and  14)  and  uncorrected 
angle  of  attack  (au)  are  plotted  on  the  lift  and  drag  curves  of  Figures  40  and  41. 

From  the  data  on  Figures  40  and  41,  it  is  seen  that  both  the  correction  procedures 
tend  to  decrease  the  correlation  with  the  so-called  cor  recti  on -free  30-  x  60-foot 
tunnel  data.  Therefore,  from  this  information,  it  would  be  hard  to  justify  the  need 
for  any  wall  corrections.  It  is  unfortunate  that  there  are  not  more  correction-free 
data  available  to  substantiate  this  finding. 

Another  interesting  result  of  Figures  40  and  41  is  that  as  the  uncorrected  lift  co¬ 
efficient  (Cj^)  increases,  the  standard  corrections  become  significantly  larger 
than  those  of  Heyson's  method.  This  is  primarily  due  to  the  large  decrease  in 
corrected  dynamic  pressure  (qc)  with  increasing  ,  which  occurs  in  Heyson's 
method.  The  reduction  in  qc  caused  an  increase  inHmth  Cl  and  Cp  ,  which  in 
effect  decreases  the  magnitude  of  Heyson's  corrections.  For  a  true  comparison  of 
the  relative  magnitudes  of  the  two  corrections,  some  sort  of  q  correction  should 
be  incorporated  in  the  standard  procedure. 

Since  the  above  correlation  was  performed  on  a  configuration  which  has  a  fairly 
small  ratio  of  model  to  tunnel  area  =  0.20),  it  was  considered  important 

to  check  this  discrepancy  with  that  obtained  on  a  configuration  having  a  larger 
value  of  Am/At.  Therefore,  the  two  correction  procedures  were  applied  to  the 
data  presented  in  Reference  54  on  a  0.6  scale  model  of  the  Biequet  941  aircraft 
tested  in  the  NASA  Ames  40-  x  80-foot  tunnel.  The  model  to  tinnel  area  ratio  of 
this  configuration  is  0.  63  (medium  span  wing). 

The  data  presented  in  Reference  57  had  already  been  corrected  by  the  standard 
procedure.  Therefore,  it  was  necessary  to  remove  the  standard  corrections  be¬ 
fore  applying  Heyson's  corrections.  The  final  corrected  lift  and  drag  data  from 
the  two  methods  are  presented  in  Figures  42  through  45.  Data  are  presented  at 
two  values  of  corrected  thrust  coefficient  and  for  two  different  combinations  of 
trailing  edge  flap  deflection. 

The  standard  wall  corrections  become  significantly  higher  than  Heyson's  as  the 
uncorrected  lift  coefficient  is  increased.  It  is  therefore  assumed  that  this  situ¬ 
ation  will  exist  regardless  of  the  ratio  of  model  area  to  tunnel  area. 
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Figure  40.  The  Effect  of  Wind  Tunnel  Wall  Corrections 
to  Lift  Coefficient  Data  From  Reference  56. 


Figure  41.  The  Effect  of  W  nd  Tunnel  Wall  Corrections 
to  Drag  Coefficient  Data  From  Reference  56. 
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Figure  42.  The  Effect  of  Wind  Tunnel  Wall  Corrections 
to  Lift  Coefficient  Data  From  Reference  57. 
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Figure  44.  The  Effect  of  Wind  Tunnel  Wall  Corrections 
to  Lift  Coefficient  Data  From  Reference  57. 
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Figure  45.  The  Effect  of  Wind  Tunnel  Wall  Corrections 
to  Drag  Coefficient  Data  From  Reference  57 
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9.0  CORRELATION  OF  THEORY  WITH  TEST  DATA 


In  order  to  determine  the  general  accuracy  of  the  aerodynamic  estimation  pro¬ 
cedures  for  lift,  longitudinal  force  and  pitching  moment  coefficients,  as  described 
in  Sections  5.0  and  6.0  of  this  report,  comprehensive  data  correlations  have  been 
carried  out  and  are  presented  in  Reference  58.  The  purposes  of  this  section  are  to 
present  selected  correlation  curves  and  to  provide  explanations  for  the  major  dis¬ 
crepancies.  Several  types  of  high-lift  configurations  were  examined,  which  in¬ 
cluded  plain,  single-slotted  and  double-slotted  trailing  edge  flaps,  and  leading 
edge  slats.  Selected  results  obtained  on  these  configurations  are  presented  for 
various  thrust  coefficients  and  for  wing  angles  of  attack  up  to  50  degrees. 


Before  proceeding  with  the  discussion  of  the  correlation  curves,  it  is  important  to 
understand  the  limitations  inherent  in  the  estimation  procedures  because  of  the 
basic  two-dimensional  data  which  form  the  foundation  of  the  procedures.  This 
two-dimensional  data  base  was  obtained  from  selected  two-dimensional  model 
tests,  and  therefore  represents  only  an  average  of  the  vast  amount  of  data  avail¬ 
able.  It  is  therefore  logical  to  assume  that  the  estimation  procedures  will  provide 
correlation  only  to  the  same  degree  of  accuracy  found  in  the  basic  two-dirnensional 
data,  and  will  not  correlate  with  every  set  of  data  examined.  The  important  thing 
to  note  is  that  the  two-dimensional  data  base  has  been  consistently  derived  for  all 
types  of  high-lift  systems,  and  should  thus  provide  a  satisfactory  basis  for  selecting 
one  high-lift  system  over  another. 


Correlations  of  a  single-slotted  trailing  edge  flap  configuration  are  presented  in 
Figures  46  through  54  for  flap  deflections  of  0,  20,  and  50  degrees.  Thrust  co¬ 
efficients  of  0,  0.6,  0.9,  and  1.0  are  shown  for  each  of  these  flap  deflections. 

With  the  exception  of  the  pitching  moment  curves  for  the  20-  and  50-degree  flap 
deflections,  all  of  the  correlations  are  quite  good.  A  detailed  analysis  of  this 
pitching  moment  discrepancy  indicated  that  the  two-dimensional  increment  in  Cm 

nio 

due  to  flap  deflection,  as  determined  from  Section  2.0,  appears  to  be  too  large. 
Comparison  of  the  predicted  two-dimensional  Cm  increment  with  the  experimental 
data  of  References  20  and  28  for  large  chord  traifing  edge  flaps  shows  the  predicted 
level  to  be  correct.  However,  a  similar  comparison  with  the  two-dimensional  data 
of  Reference  2  shows  poor  agreement  for  flap  chord  to  wing  chord  ratios  greater 
than  35  percent.  Therefore,  it  is  possible  that  the  aerodynamic  estimation  pro¬ 
cedure  should  have  factors  for  large  values  of  flap  chord  to  wing  chord  ratio  in 
order  to  provide  the  necessary  correlation  tolerance. 


The  double-slotted  trailing  edge  flap  correlations  are  presented  in  Figures  55 
through  57  for  a  flap  deflection  of  40  degrees  at  thrust  coefficients  of  0. 6,  0. 8, 
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and  0,9.  The  significant  discrepancy  in  this  data  again  lies  in  the  area  of  Cm 
prediction.  In  this  case,  the  two-dimensional  estimate  is  conservative  rathei^ 
than  overpredicted.  This  further  points  out  the  strong  need  of  additional  work  in 
this  area  of  the  moment  estimation  procedure. 


The  correlation  data  for  the  leading  edge  slat  configuration  in  conjunction  with 
the  50-degree  single-slotted  trailing  edge  flap  are  shown  in  Figures  58  through  60. 
Once  again,  thrust  coefficients  of  0,  0.6,  0.9,  and  1 . 0  are  presented.  In  addition 
to  the  pitching  moment  discrepancy,  which  was  discussed  in  the  previous  para¬ 
graph,  there  is  a  large  discrepancy  in  the  lift  curves  at  the  low  angles  of  attack. 
This  discrepancy  results  from  the  typical  problem  of  leading  edge  devices:  flow 
separation  at  low  angles  of  attack.  The  estimation  procedure,  as  it  now  exists, 
does  not  treat  the  separated  flow  phenomena  whicivare  predominant  on  this  type  of 
configuration.  According  to  the  test  data,  the  addition  of  this  leading  edge  slat 
actually  decreases  the  available  maximum  lift  coefficient  of  the  total  system.  This 


is  the  primary  reason  for  the  rather  poor  C 


lMAX 


correlation  on  this  configuration. 
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Figure  46.  Total  Lift  Coefficient  Versus  Wing  Angle  of  Attack,  6^  =  0° . 
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Figure  47.  Total  Lift  Coefficient  Versus  Total  Drag  Coefficient, 
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Figure  55.  Total  Lift  Coefficient  Versus  Wing  Angle  of  Attack,  6^  =  40°. 
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Figure  56.  Total  Lift  Coefficient  Versus  Total  Drag  Coefficient, 
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Figure  58.  Total  Lift  Coefficient  Versus  Wing  Angle  of  Attack, 
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Figure  59.  Total  Lift  Coefficient  Versus  Total  Drag  Coefficient, 


10.0  PERFORMANCE  PROGRAM  TO  CALCULATE  TAKEOFF, 

LANDING  AND  TRANSITION  OF  V/STOL  AIRCRAFT 

The  digital  computer  program  described  herein  was  developed  to  evaluate  the  effect 
of  various  high-lift  systems  on  takeoff,  landing  and  transition  performance.  Com¬ 
plete  sets  of  program  listings,  program  flow  charts,  and  input  data  for  a  sample 
case  are  included  in  Reference  61.  The  summary  equations  in  Sections  4.0  and 
5.0  have  been  programmed  (Reference  45)  to  supply  aerodynamic  data  input  to  the 
performance  program. 

iO.  1  EQUATIONS  OF  MOTION 

The  equations  used  in  this  calculation  system  are  derived  for  an  earth  axis  motion 
system  (ground  speed,  vertical  speed).  Acceleration  in  the  vertical  speed  direc¬ 
tion  is  given  by 

Az  =  (g/W)  jlift  cosy  +  Drag  siny-W  -  (EQ)T  sin  (C*w  +  iT  +  y)j  (158) 

and  acceleration  in  the  ground  speed  direction  is  given  by 

Ax  =  (gAV)  j^Drag  cosy- Lift  sin y-(EQ)Tcos  (aw  +  iT  +  y)J+Az)i  (159) 

where  AX  =  horizontal  acceleration 
Az  =  vertical  acceleration 

Drag  =  summation  of  drag  and  thrust  forces  along  the  free-stream 
velocity  vector  (positive  forward) 

EQ  =  flight  condition  index 

=  -1.0  for  hovering  or  descending  flight 
=  +1.0  for  thrust  reversal 
=  0  for  all  other  conditions 

g  =  acceleration  due  to  gravity 

i^,  =  thrust  incidence 

Lift  =  summation  of  lift  and  thrust  forces  perpendicular  to  the  free- 
stream  velocity  vector  (positive  up) 

T  =  thrust 

W  =  weight 

12  i 


0<w  =  wing  angle  of  attack 

y  =  flight  path  angle 

M  =  coefficient  of  friction  (ground  roll  portions  only) 
The  forces  and  angles  are  shown  in  Figure  61. 


Figure  61.  Force  and  Angle  Diagram  for  Equations  of  Motion. 


10.2  PROGRAM  DESCRIPTION 

The  primary  objectives  during  the  development  of  this  calculation  system  were: 

1 .  Minimum  input  data  generation 

2.  Reasonable  piloting  techniques 

3.  Ability  to  refine  the  subroutines  with  minimum  interfacing  problems 

These  objectives  were  satisfied  by  designing  a  modular  program  with  the  five 
following  functional  groups: 

1.  The  main  program  (direct)  is  the  executive  routine  and  controls  the 
reading  of  inputs,  initializing  weights,  altitudes,  and  temperatures,  and 
the  executing  of  the  configuration  subprograms  for  a  given  weight,  alti¬ 
tude  and  temperature  matrix. 

2.  The  configuration  subprograms  are  designed  to  "drive"  and  monitor  the 
particular  calculation  desired,  e.g.,  deflected-slipstream  short  takeoff, 
tilt -wing  vertical  landing. 

3.  The  physical  data  subroutines  are  used  to  store  and  retrieve  aerodynamic 
date,  propulsion  characteristics,  and  atmospheric  properties.  In  addi¬ 
tion,  an  aerodynamic  subroutine  calculates  accelerations  acting  on  the 
aircraft. 

4.  The  general  use  subroutines  provide  calculations  or  timing  cues  used  in 
one  or  more  of  the  configuration  subroutines. 

5.  The  support  subroutines  do  table  look-ups  and  evaluate  curve  fits. 

The  program  can  perform  calculations  for  deflected-slipstream  takeoffs  and  land¬ 
ings  and  tilt-wing  short  takeoffs,  vertical  takeoffs,  vertical  landings,  and  transi¬ 
tions.  Tilt-prop  V/STOL  and  tilt-wing  short  landing  calculation  subprograms 
have  not  been  included  in  this  system.  Provisions  have  been  made,  however,  in 
the  physical  data  and  general  use  subroutines  to  support  calculations  for  any 
propeller-driven  V/STOL  aircraft. 

In  addition  to  the  calculation  system,  a  trimming  program  was  written.  A  sche¬ 
matic  of  the  program  flow  and  structure  is  given  in  Figure  62. 
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wructure  by  Functional  Group. 


10.3  SIMPLIFYING  ASSUMPTIONS 


The  major  assumptions  made  to  simplify  the  calculations  and  trajectories  are 
listed  by  subroutine  or  program  name.  The  assumptions  reflect  reasonable  pilot¬ 
ing  technique  and  reduce  a  three-degree-of -freedom  problem  to  a  more  manage¬ 
able  two-degree -of -freedom  problem.  In  all  cases  the  assumptions  can  be  elim¬ 
inated  by  substituting  a  more  refined  subroutine  for  the  routine  presently  in  the 
system. 

10.3.1  Main  Program 

The  assumptions  in  the  main  program  are: 

1.  The  tail-off  pitching  moment  is  trimmed  using  the  expression 

CmTAiL-OFF+ ClTAlL-OFF(°* 25_CG) 
CLtRIMMED_ClT  AIL-OFF  +  i 

H/c 


where  4  is  the  nondimensional  tail  arm. 
H/c 


CG  is  the  location  of  the  center  of  gravity  in  relation  to  the  mean 
aerodynamic  chord. 

CjYirp  ATT  is  referenced  to  the  0. 25  c  station. 

I  AIL-OFF 


2.  Incremental  drag  due  to  trimming  is  several  orders  of  magnitude  smaller 
than  the  total  aircraft  drag  and  is  considered  to  be  negligible. 

3.  All  trimming  on  the  tilt-wing  configuration  is  with  the  tail  rotor. 

4.  Atmospheric  properties  are  not  recomputed  during  the  problem.  Due 
to  the  small  altitude  changes  encountered,  runway  altitude  and  ambient 
temperature  are  used  in  all  density  expressions. 

5.  Aircraft  gross  weight  is  not  reduced  as  fuel  is  used.  The  time  involved 
in  any  of  the  maneuvers  is  short,  and  even  at  maximum  power,  fuel  used 
is  less  than  0.5%  of  the  aircraft  gross  weight. 

6.  Zero  wind  conditions  are  assumed  by  all  routines  in  this  system. 
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10.3.2  Configuration  Subprograms 

The  assumptions  in  the  configuration  subprograms  are: 

Deflected-Slipstream  Short  Landing  —  The  trajectory  assumed  for  this  configuration 
•  is  represented  in  Figure  63. 

1.  "Start  at  the  Obstacle"  —  Airspeed  at  this  point  is  a  percentage  of  power- 
on  stall  speed.  Sink  speed  is  an  input.  Both  airspeed  and  sink  speed  are 
held  constant  to  touchdown. 

2.  ’Touchdown"  —  If  the  aircraft  attitude  at  touchdown  will  cause  a  "nosewheel- 
first"  landing,  the  aircraft  is  instantaneously  rotated  to  ground  attitude  at 
touchdown. 

3.  "Braking  Devices  On"  —  Braking  and  rolling  coefficients  of  friction  are 
not  a  function  of  velocity. 

Deflected-Slipstream  Short  Takeoff  —  The  trajectory  assumed  for  this  configura¬ 
tion  is  represented  in  Figure  64. 

1.  "Start  Takeoff"  —  The  entire  maneuver  is  performed  at  a  constant  power 
setting.  Rolling  coefficient  of  friction  is  not  a  function  of  velocity. 

2.  "Rotation  and  Lift-Off"  —  Rotation  speed  is  a  percentage  of  power-on 
stall  speed.  The  aircraft  is  rotated  at  a  constant  angular  rate  until  the 
maximum  angle  of  attack  is  attained. 

3.  "Clear  the  Obstacle"  —  If  the  maximum  angle  of  attack  is  attained  before 
the  obstacle  height  is  reached,  the  maximum  angle  is  held  over  the 
obstacle. 

Tilt- Wing  Short  Takeoff  —  The  trajectory  assumed  for  this  configuration  is  rep¬ 
resented  in  Figure  65. 

1.  "Start  Takeoff"  —  Wing  tilt  angle  is  set  at  the  beginning  of  the  maneuver. 

Flap  deflection  is  a  sine  function  of  wing  tilt  angle  (e.g. ,  wing  tilt  =  0, 
flap  =  0;  wing  tilt  =  45  deg,  flap  is  at  the  maximum  deflection;  wing  tilt  = 

90  deg,  flap  =  0).  Coefficient  of  rolling  friction  is  not  a  function  of 
velocity. 

2.  "Airborne"  —  Fuselage  angle  of  attack  and  wing  tilt  angle  are  constant 
throughout  the  maneuver.  Tail  rotor  thrust  is  in  the  lift  direction  for 
all  wing  tilt  angles.  (This  is  a  predetermined  percentage  of  total  power 
as  a  function  of  tilt  angle  and  is  used  to  replace  the  trim  equation  of 
10.3.1). 


3.  "Cleat  the  Obstacle"  —  The  problem  is  terminated  at  the  obstacle.  Once 
the  aircraft  is  above  the  obstacle,  forward  speed  is  sufficient  for  an  easy 
transition  to  conventional  flight. 

Tilt-Wing  Vertical  Landing  —  The  vertical  landing  trajectory  is  shown  in  Figure  66. 
To  attempt  a  vertical  landing,  the  aircraft  must  have  a  thrust  to  weight  ratio  greater 
than  1.07  to  insure  an  adequate  maneuver  margin. 

1.  "Conventional  Flight"  —  The  maneuver  is  started  at  a  steady-state  speed 
which  is  a  percentage  of  power-off  stall  speed. 

2.  'Transition"  —  The  wing  is  tilted  at  a  constant  angular  rate.  Flap  deflec¬ 
tion  is  a  sme  function  of  wing  tilt  angle.  Zero  flight  path  angle  is  main¬ 
tained  throughout  the  transition.  Tail  rotor  thrust  is  in  the  lift  direction 
for  ail  wing  tilt  angles. 

3.  "Hover"  —  Maximum  wing  angle  of  attack  during  this  maneuver  is  90 
degrees  plus  the  angle  of  zero  lift. 

4.  'Touchdown"  —  Input  sink  speed  is  maintained  with  throttle  control  during 
the  vertical  descent. 

Tilt-Wing  Vertical  Takeoff  —  The  vertical  takeoff  trajectory  is  shown  in  Figure  67. 
To  attempt  a  vertical  takeoff,  the  aircraft  must  have  a  thrust  to  weight  ratio  greater 
than  1. 05  to  insure  an  adequate  maneuver  margin. 


1.  "Start  Takeoff"  —  Takeoff  is  initiated  at  thrust  equals  weight.  The  throttle 
is  then  advanced  to  maximum  thrust. 

2.  "Vertical  Ascent"  —  Fuselage  pitch  angle  is  held  constant  during  this 
phase.  Tail  rotor  thrust  is  in  the  lift  direction  for  all  wing  tilt  angles. 
This  phase  terminates  when  a  minimum  rate  of  climb  and  a  minimum 
altitude  are  exceeded. 

3.  "Start  Transition  -  Pitch  Fuselage  Down"  —  The  fuselage  nose  is  pitched 
down  to  a  preset  angle  to  make  the  transition  to  conventional  flight 
quicker. 

4.  "Transition  -  Tilt  Wing  Down"  —  The  wing  is  tilted  down  at  a  constant 
angular  rate.  Flap  deflection  is  a  sine  function  of  wing  tilt  angle. 

5.  "Conventional  Flight"  —  The  problem  terminates  when  wing  tilt  angle 
reaches  zero. 
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10.3.3  Physical  Data  Subroutines 

The  ass  -mptione  in  the  physical  data  subroutines  are: 

Aerodynamic  Data  Subroutine 

1.  When  the  aircraft  is  in  a  hovering  or  vertical  descent  mode  (CT  =  1.00), 
Cl  and  Cp  equal  zero. 

2.  Aerodynamic  coefficients  include  both  indirect  and  direct  thrust  effects. 

3.  For  reverse  thrust,  wing  lift  equals  zero. 

Atmospheric  Properties  Subroutine 

The  standard  day  temperature  lapse  rate  is  valid  for  all  temperature/altitude 

combinations  to  be  considered. 

Propulsion  Characteristics  Subroutine 

1.  One  propeller  per  engine. 

2.  A  constant  percentage  of  optimum  propeller  efficiency  is  valid  for  the 
speed  regime  considered. 

3.  For  tilt-wing  aircraft,  tail  rotor  power  extraction  is  a  function  of  wing 
tilt  angle  only. 

4.  Input  data  are  valid  throughout  the  altit  .de  change  during  the  maneuver. 

5.  Axial  inflow  to  the  propeller.  A  sample  calculation  was  made  with  oblique 
inflow  and  compared  with  the  axial  inflow  case.  The  difference  was 
negligible. 

10.3.4  General  Use  Subroutines 

The  assumptions  in  the  general  use  subroutines  are: 

Subroutine  BRAKE  (provides  timing  cues  for  ground  roll  portion  of  short  landings) 

1.  Aircraft  nose  is  rotated  to  a  three-point  attitude  at  a  constant  rate. 

2.  Brake  application  occurs  1  second  after  touchdown. 

3.  To  use  reverse  thrust,  the  engines  are  throttled  back  to  idle  power  before 
reverse  thrust  is  initiated. 

4.  The  wing  is  tilted  down  at  a  constant  rate  after  touchdown  for  all  tilt-wing 
and  tilt-prop  aircraft. 
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Integration  Subroutine 

Trapezoidal  rule  integration  with  a  small  time  step  (approximately  0.25  sec) 
is  sufficiently  accurate. 


Power-On  Stall  Speed  Convergence 


C 


lmax 


occurs  at  an  angle  of  attack  determined  by  pilot  comfort  rather  than 


aerodynamic  separation. 
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11.0 


CONCLUSIONS 


1.  The  procedures  to  predict  the  two-dimensional  aerodynamic  characteristics  of 
high-lift  devices  furnish  reasonable  estimates,  generally  within  10  percent  of 
experimental  results,  that  can  be  used  in  the  span  load  program.  The  trend 
from  General  Dynamics  experimental  data  indicates  that  the  maximum  lift 
levels  from  NACA  data  are  approximately  10  to  20  percent  lower.  This  may 
be  attributed  to  the  General  Dynamics  testing  technique  which  utilizes  side- 
wall  blowing  to  eliminate  adverse  wall  boundary  layer  interference  effects. 

2.  The  program  initiated  at  Convair  to  develop  a  more  basic  approach  to  the  prob¬ 
lem  of  estimating  the  lift  and  longitudinal  force  coefficients  was  well  justified. 
It  is  considerably  less  limited  than  existing  empirical  methods  in  that  it  re¬ 
quires  only  two-dimensional  data  for  application  to  general  configurations. 

The  procedures  for  predicting  lift  and  longitudinal  force  coefficients  of  a  wing- 
flap  combination  give  satisfactory  results  at  all  thrust  coefficients.  However, 
the  procedures  for  predicting  pitching  moment  coefficients  result  in  erratic 
correlations  for  flaps -deflected  cases. 

3.  Tne  entire  task  of  correlating  wind  tunnel  data  that  have  been  corrected  for 
wall  effects  with  actual  flight  test  conditions  could  not  be  accomplished  due 
to  the  lack  of  a  good  set  of  flight  test  data  for  comparison.  The  evaluation  of 
tunnel  wall  effects  indicated  that  currently  available  correction  procedures  for 
lift  and  drag  yield  erroneous  results  and  should  not  be  applied. 

4.  The  complexity  of  the  performance  program  made  numerous  assumptions  and 
capability  restrictions  desirable.  For  example,  the  equations  of  motto?'  were 
reduced  to  two  dimensions  (i.e.,  no  pitch  dynamics),  and  several  terms,  such 
as  pitch  rates,  wing  tilt  angular  rates,  and  pilot  braking  response  time,  were 
used  as  predetermined  inputs.  These  assumptions  were  made  to  expedite  the 
development  of  the  system;  however,  the  accuracy  of  the  methods  employed 

is  sufficient  to  reflect  the  impact  of  lugh-lift  devices  on  V/STOL  aircraft. 

The  program  is  a  viable  tool  for  performance  estimation. 
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12.0  RECOMMENDATIONS 


1.  The  problems  associated  with  predicting  the  aerodynamic  forces  acting  on  a 
wing  immersed  in.  a  slipstream  should  be  investigated  further.  To  improve 
the  basic  procedures,  it  is  necessary  that  (1)  additional  correlations  of  pitch¬ 
ing  moments  be  made  with  experimental  data  to  develop  improved  empirical 
factors  for  large  chord  ratio  flaps,  (2)  methods  be  developed  to  incorporate 
the  pronounced  nonlinear  characteristics  of  leading  edge  devices  at  low  angles 
of  attack,  (3)  the  downwash  estimating  procedure  be  programmed  for  digital 
computer  application  and  correlations  be  made  with  experimental  data,  and 
(4)  tilt-prop  configuration  methodology  be  developed  and  incorporated  into  the 
aerodynamic  and  performance  programs. 

2.  Jet  flap  and  boundary  layer  control  devices  should  be  investigated  so  that 
methods  can  be  developed  to  describe  the  effects  of  these  devices  on  a  wing 
immersed  in  a  slipstream. 

3.  A  study  should  be  performed  to  define  a  wind  tunnel/flight  test  program  to 
correlate  V/STOL  vehicle  aerodynamic  characteristics . 
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APPENDIX 


ANALYSIS  OF  THE  EFFECTS  OF  HIGH-LIFT  DEVICES  ON 
THE  PERFORMANCE  OF  A  TILT-WING  CONFIGURATION 

The  analysis  of  the  effects  of  high-lift  devices  on  the  performance  of  a  tilt-wing 
configuration  similar  to  the  CL-84  is  discussed  in  this  section.  The  general 
arrangement  of  the  aircraft  (powered  by  two  T64-GE-7  engines)  is  shown  in  Figure 
68,  The  geometry  is  listed  below: 


Wing 

Horizontal  Tail 

s 

=  291  sq  ft 

S  =  95  sq  ft 

AR 

=  5.18 

t/c  =0.14 

b 

-  38.8  ft 

c 

=  7.5  ft 

\ 

=  1.0 

Airfoil  NACA  63g  -  418 

» ,  i 

Fuselage 

Vertical  Tail 

=  44  ft 

Center  Outboard 

db 

SWET 

=  7.7  ft 

S  =  33  sq  ft  S  =27  sq  ft 

=  764  sq  ft 

t/c  =0.14  t/c  =0.12 

Nacelle 

Propeller 

£n 

=  14.24  ft 

Dp  =16.5  ft 

SB 

=  13  sq  ft 

Number  of  Blades  =  4 

The  methods  discussed  in  Section  2.0  were  used  to  estimate  the  two-dimensional 
characteristics  of  the  35-percent-chord  trailing  edge  flaps  shown  in  Figure  69. 
These  characteristics  plus  wing  geometry  were  inputs  into  the  span  load  program. 
The  results  from  the  span  load  program  plus  the  nacelle  drag  were  then  used  in 
the  summary  equations  of  Sections  5.0  and  6,0  to  generate  the  tail -off  aerodyna¬ 
mic  characteristics.  These  data,  fuselage  drag,  vertical  tail  drag,  horizontal  tail 


144 


drag,  geometry,  and  weights,  were  inputs  into  the  performance  program  described 
in  Section  10,0.  The  performance  for  the  aircraft  with  the  various  high-lift  de¬ 
vices  is  shown  in  Figures  70  through  75, 

The  flight  paths  for  vertical  takeoff  to  conversion  (wing  tilt  angle  at  4.16  degrees) 
for  a  gross  weight  of  20,000  pounds  are  shown  in  Figure  70.  The  35-percent- 
chord  trailing  edge  flaps  furnished  sufficient  lift  to  allow  the  vehicle  to  perform 
the  transition  from  vertical  takeoff  to  conversion  up  to  the  maximum  vertical  take¬ 
off  weight.  The  selection  of  a  trailing  edge  flap  has  a  pronounced  effect  on  ground 
distance  and  the  altitude  at  which  the  conversion  is  completed.  Typical  distances, 
altitudes,  and  velocities  to  conversion  for  a  range  of  gross  weights  up  to  the  maxi¬ 
mum  vertical  takeoff  gross  weight  of  25,000  pounds 

are  shown  in  Figure  71.  The  overall  effect  of  improving  the  lifting  capability  of 
the  trailing  edge  flap  is  to  drastically  reduce  the  distance  to  conversion  and  allow 
the  aircraft  to  maintain  an  altitude  at  conversion  which  is  almost  constant  with 
weight. 

The  effects  of  wing  tilt  angle  on  the  short-field  takeoff  distances  at  a  gross  weight 
of  24,000  pounds  are  shown  in  Figure  72.  Wing  tilt  angles  of  40  degrees  for  the 
double -slotted  flap  and  45  degrees  for  the  plain  and  single -slotted  flap  configura¬ 
tions  gave  the  shortest  distances  over  the  obstacle.  These  wing  tilt  angles  were 
used  for  the  short-field  takeoff  analysis . 

The  effects  of  trailing  edge  flaps  on  the  short-field  takeoff  are  shown  in  Figure  73. 
The  wing  is  tilted  to  a  prescribed  angle  and  the  trailing  edge  flaps  are  fully  de¬ 
flected  to  keep  the  ground  distances  down  to  100  feet.  The  results  for  the  trailing 
edge  flaps  analyzed  are  pronounced.  For  a  500-foot  total  distance  requirement  for 
an  overload  condition,  the  single-  and  double -slotted  flap  configurations  would 
allow  an  additional  2200  and  4100  pounds  of  payload  respectively. 

The  double -slotted  flap  configuration  was  then  analyzed  for  the  vertical  landing  and 
deflected-slipstream  takeoff  performance  to  demonstrate  the  capability  of  the  pro¬ 
gram.  Time  and  distance  for  a  typical  vertical  landing  for  the  double-slotted  flap 
are  shown  in  Figure  74.  The  maximum  vertical  landing  weight  is  24,550  pounds. 
The  major  effect  of  reducing  the  gross  weight  is  a  reduction  in  time  and  total  dis¬ 
tance  to  hover.  Typical  time  required  for  tire  maneuvers  is  10  seconds. 

Another  20  seconds  is  required  from  hover  to  touchdown. 

Deflected -slipstream  takeoff  performance  for  the  double -slotted  flap  is  shown  in 
Figure  75.  The  wing  is  not  tilted  but  the  flap  is  fully  deflected  for  this  maneuver. 
The  distances  at  the  lower  gross  weights  appear  low  when  compared  to  the  Brequet 


max.  available  thrust  \ 
lT05  / 
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and  Charger  takeoff  data.  The  following  table  shows  that  the  tilt -wing  configura¬ 
tion  selected  for  analysis  has  substantially  higher  thrust  to  weight  ratio  and  wing 
loading  than  the  Brequet  or  Charger.  The  analysis  has  served  as  a  very  satis¬ 
factory  check  of  the  performance  program.  A  further  check  should  be  made  with 
the  Brequet  or  the  Charger  to  show  the  full  program  capability. 


COMPARISON  OF  DEFLECTED-SLIPSTREAM  PERFORMANCE 

AIRCRAFT 

TILT-WING 

BREQUET 

CHARGER 

Weight  (lb) 

29,100 

48,500 

7,800 

S  (sq  ft) 

291 

902 

191 

W/S 

100 

54 

41 

No.  of  Engines 

2 

4 

2 

SHP—.—. _ 

TOTAL 

6,870 

6,000 

1,300 

W/SHP___a  _ 

TOTAL 

4.23 

8.08 

6.0 

Ground  Distance 

240 

610 

275 

Total  Distance  Over 

50-Ft  Obstacle 

640 

985 

610 

Gross  Weight  -  20,000  lb 


/ 


Gingle-Slotted 

Flap 


iround  Distance  -  ft 


Figure  70.  Typical  Vertical  Takeoff  Flight  Paths. 
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Gross  Weight  -  1000  lb 

Figure  73.  Short-Field  Takeoff  Distances  With  Wing  Tilted  Versus  Gross 


Total  Horizontal  Distance 
to  Vertical  Landing  -  ft 


Gross  Weight  -  1000  lb 


Figure  74.  Distance  and  Time  for  Vertical  Landing  Versus  Gross  Weight. 
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